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Summary 


In  this  paper  explicit  necessary  and  sufficient  conditions  are  estab¬ 
lished  for  the  ordinary  and  strong  ellipticity  of  the  three-dimensional 
field  equations  in  the  nonlinear  equilibrium  theory  of  incompressible, 
homogeneous  and  isotropic,  hyperelastic  solids.  The  resulting  system  of 
inequalities  involves  the  local  principal  stretches  directly  and  in  addition 
restricts  the  first  and  second  partial  derivatives  of  the  strain-energy 
density  with  respect  to  the  deformation  invariants  or  the  principal  stretches 
The  conditions  of  ordinary  and  strong  ellipticity  are  found  to  coalesce  for 
materials  that  obey  the  Baker-Ericksen  inequalities  and  possess  a  positive 
shear  modulus  at  infinitesimal  deformations.  Various  implications  of  these 
ellipticity  conditions  for  special  classes  of  materials  and  deformations  are 
explored. 


The  results  communicated  in  this  paper  were  obtained  in  the  course  of  an 
investigation  supported  in  part  by  Contract  N00014-75-C-0196  with  the  Office 
of  Naval  Research  in  Washington,  D.C. 


Introduction 

Issues  related  to  the  ellipticity  of  the  equations  governing  the  fi¬ 
nite  equilibrium  theory  of  perfectly  elastic  solids  and  to  the  possible 
breakdown  of  ellipticity  in  nonlinear  elastostatics,  have  attracted  in¬ 
creasing  attention  during  recent  years.  Moreover,  this  interest  has 
sprung  from  diverse  and  -  to  some  extent  -  conflicting  motivations. 

Thus,  some  of  the  work  to  which  we  are  alluding  seeks  to  extend  the 
scope  of  the  theory  to  equilibrium  solutions  of  reduced  regularity  that 
encompass  discontinuous  deformation  gradients  of  the  kind  associated  with 
so-called  localized  shear  failures.  The  emergence  of  such  singular  equi¬ 
librium  fields  in  homogeneous  hyperelastic  solids  is  accompanied 
by  a  failure  of  ordinary  ellipticity.  Closely  allied  in  purpose  are 
investigations  concerned  with  bifurcations  of  equilibrium  solutions  that 
are  contingent  upon  a  loss  of  strong  ellipticity.  In  contrast,  other 
related  work  has  a  different  incentive:  it  aims  chiefly  at  restrictions 
of  the  strain-energy  density  arising  from  the  postulate  of  "material  sta¬ 
bility"  in  the  sense  of  Hadamard,  which  precludes  a  loss  of  strong  ellip¬ 
ticity  in  the  elastostatic  field  equations. 

Although  the  present  paper  bears  on  both  of  these  objectives,  it 
derives  its  impetus  from  the  first  of  the  foregoing  two  motivations  and 
continues  a  sequence  of  studies  initiated  in  [1],  The  latter  was  prompted 
by  the  surmise  that  a  certain  crack  problem  in  finite  elastostatics  of 
compressible  hyperelastic  solids  fails  to  admit  a  solution  of  unlimited 
smoothness  for  a  particular  hypothetical  isotropic  material.  This  con¬ 
jecture,  in  turn,  suggested  that  the  corresponding  displacement  equations 
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of  equilibrium  suffer  a  loss  of  ordinary  ellipticity  in  the  presence  of 
severe  enough  deformations.  That  such  Is  indeed  the  case  is  borne  out 
by  the  analysis  in  [1],  which  supplies  an  appropriate  necessary  and  suf¬ 
ficient  criterion  of  ellipticity. 

Explicit  conditions  necessary  and  sufficient  for  ordinary  and  strong 
ellipticity  of  the  two-dimensional  field  equations  in  the  theory  of  finite 
plane  strain  for  compressible  isotropic  hyperelastic  bodies  are  deduced 
in  [2].  The  results  of  [2]  are  applied  in  [3]  to  a  local  study  of  plane 
deformation  fields  that  possess  continuous  displacements,  but  exhibit 
finite  jump  discontinuities  in  the  deformation  gradient.  The  energetics 
of  such  "elastostatic  shocks"  are  further  explored  by  Knowles  [4]. 

A  necessary  and  sufficient  ellipticity  criterion  confined  to  anti- 
plane  shear  deformations  of  a  class  of  incompressible  isotropic,  perfectly 
elastic  materials  is  included  in  [4]  and  applied  in  [5].  Further,  results 
essentially  analogous  to  those  contained  in  [2],  [3],  [4],  but  pertaining 
to  plane  deformations  of  incompressible  hyperelastic  solids,  are  derived 
by  Abeyaratne  [6],  who  limits  his  attention  to  ordinary  ellipticity. 

Specific  boundary-value  problems  involving  a  loss  of  ellipticity  and 
the  concomitant  appearance  of  elastostatic  shocks  are  treated  in  a  number 
of  publications,  all  but  one  of  which  deal  asymptotically  with  the  equi¬ 
librium  field  near  the  tip  of  a  crack  in  an  incompressible  body  subjected 
to  anti -plane  shear.  References  to  these  papers  can  be  found  in  a  recent 
survey  [7]J 

The  present  Investigation  furnishes  explicit  necessary  and  sufficient 

^To  these  we  append  an  asymptotic  study  by  Abeyaratne  [8],  which  has  ap¬ 
peared  since. 
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elllptlcity  conditions  appropriate  to  the  general  three-dimensional  field 
equations  for  incompressible  isotropic  hyperelastic  bodies.  Further,  both 
ordinary  and  strong  ellipticity  are  considered,  and  it  is  shown  that  the 
corresponding  ellipticity  criteria  merge  for  a  material  obeying  the 
Baker-Ericksen  Inequalities,  provided  its  shear  modulus  is  positive  at 
infinitesimal  deformations.  Thus,  In  these  circumstances,  ordinary  implies 
strong  ellipticity. 

As  far  as  additional  related  literature  is  concerned,  we  cite  first 
an  analysis  of  localized  shear  failures  due  to  Rudnicki  and  Rice  [9],  as 
well  as  a  bifurcation  analysis  by  Hill  and  Hutchinson  [10].  Both  of  these 
papers  presuppose  a  potential  loss  of  ellipticity  in  equilibrium  continuum 
mechanics;  their  constitutive  settings,  however,  go  beyond  purely  elastic 
behavior.  Finally,  it  should  be  emphasized  that  the  work  reported  here 
is  rather  closely  connected  with,  and  complementary  to,  a  series  of  in¬ 
vestigations  due  to  Sawyers  and  Rlvlln,  starting  with  [11].  These  papers 
are  c1t*d  In  a  survey  article  by  Sawyers  [12],  which  summarizes  the  results 
obtained.  The  latter  include  various  necessary  conditions  for  strong 
ellipticity  within  the  context  of  the  equilibrium  theory  of  incompressible, 
homogeneous  and  isotropic,  hyperelastic  solids.  These  conditions,  which 
are  sufficient  merely  for  certain  restricted  classes  of  such  materials  or 
when  the  deformation  exhibits  a  particular  degeneracy,  are  included  among 
the  complete  ellipticity  criteria  established  in  the  present  paper. 

In  Section  1,  which  Is  partly  expository,  we  assemble  some  essential 
Ingredients  of  the  nonlinear  equilibrium  theory  for  homogeneous  incompres¬ 
sible  hyperelastic  bodies.  Here  we  also  define  the  appropriate  concepts 
of  ordinary  and  strong  ellipticity.  Proceeding  from  these  two  definitions, 
we  establish  corresponding  necessary  and  sufficient  ellipticity  conditions 


j  Atnt:  -in  _  _ 
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that  encompass  —  but  are  not  confined  to  —  the  special  case  of  material 
isotropy  treated  in  the  remainder  of  the  paper. 

Section  2  is  devoted  to  the  derivation  of  explicit  necessary  condi¬ 
tions  for  ordinary  ellipticity  pertaining  to  the  subclass  of  isotropic 
materials.  The  ensuing  inequalities  involve  the  principal  stretches  di¬ 
rectly,  as  well  as  through  the  first  and  second  partial  derivatives  of 
the  strain-energy  density  with  respect  to  the  deformation  invariants. 
Further,  at  the  end  of  Section  2,  we  deduce  an  equivalent  system  of  in¬ 
equalities  necessary  for  ordinary  ellipticity  that  depends  upon  the 
elastic  potential  exclusively  through  its  first  and  second  gradients  with 
respect  to  the  principal  stretches  and  is  fully  symmetric  in  the  latter. 

In  Section  3  we  first  establish  the  sufficiency  of  the  necessary 
conditions  for  ordinary  ellipticity  arrived  at  in  Section  2.  We  then  de¬ 
rive  explicit  necessary  and  sufficient  conditions  for  the  strong  ellip¬ 
ticity  of  the  elastostatic  field  equations  on  the  assumption  of  material 
isotropy.  Finally,  we  show  that  the  conditions  of  ordinary  and  strong 
ellipticity  coalesce  in  this  instance  under  certain  mild  and  physically 
plausible  additional  restrictions  of  the  material  response. 

The  concluding  Section  4  deals  with  applications  of  the  ellipticity 
conditions  established  earlier  to  particular  classes  of  deformations  and 
to  special  types  of  elastic  materials  within  the  category  under  consider¬ 
ation.  Here  we  examine  the  degenerate  instances  of  a  locally  axisymmetric 
and  a  locally  plane  deformation,  as  well  as  the  case  in  which  the  elastic 
potential  depends  merely  on  a  single  deformation  invariant.  Next,  we 
Infer  the  eiiipticity,  at  all  deformations,  of  the  equilibrium  field 
equations  appropriate  to  a  Mooney-Rivl in  material.  Finally,  as  an  illus¬ 
trative  example,  we  discuss  in  detail  the  domain  of  ellipticity  in  the 


space  of  the  principal  stretches  for  a  specific  material  that  has  an 
elastic  potential  of  a  form  proposed  by  Ogden  [13]. 


-6- 


1 .  Preliminaries  from  finite  elastostatics.  Ordinary  and  strong  ellip- 

ticity  for  imcompressible,  hyperelastic  solids. 

In  this  section  we  recall  certain  prerequisites  from  the  finite  equi¬ 
librium  theory  of  homogeneous  incompressible  hyperelastic  solids.  We  then 
define  the  notions  of  ordinary  and  strong  ellipticity  in  the  present  con¬ 
text  and  deduce  necessary  and  sufficient  ellipticity  conditions  that  are 
not  contingent  upon  any  material  symmetry  restrictions. 

Throughout  this  paper,  uppercase  boldface  letters  denote  second-order 
tensors  as  well  as  three-by- three  matrices;  lowercase  boldface  letters 
denote  vectors  and  also  three-rowed  column  matrices.  Further,  the  same 
boldface  letter  will  be  used  to  designate  a  tensor  or  vector  and  its  ma¬ 
trix  of  scalar  components  in  the  underlying  rectangular  Cartesian  coordi¬ 
nate  frame. 

Let  ft  be  the  three-dimensional  open  region  occupied  by  the  interior 
of  a  body  in  an  undeformed  reference  configuration.  A  deformation  of  the 
body  is  then  described  by  a  transformation 

y  =  y(x)  =  x  +  u(x)  on  ft  ,  (1.1) 

which  maps  ft  onto  a  domain  ft*.  Here  x  is  the  position  vector  of  a 
generic  point  in  ft,  y(x)  is  its  deformation  image  in  ft*,  and  u  is 
the  displacement  vector  field.  Thus,  x^  and  y^  are  the  Cartesian  ma¬ 
terial  and  spatial  coordinates,  respectively.^  We  shall  suppose  for  the 
time  being  that  the  mapping  y  is  twice  continuously  differentiable  and 


^ Latin  subscripts  have  the  range  (1,2,3)  and  summation  over  repeated  sub¬ 
scripts  is  taken  for  granted. 


one-to-one  on  ft.  Next,  let 


F  =  vy=  [3y /ax  ]  =  [6  +3u./axJ,  J  =  det  F  on  ft,  (1.2 

~  ~  i  j  ij  i  j  ~ 

so  that  £  is  the  deformation-gradient  tensor  and  J  the  Jacobian  deter¬ 
minant  (volume  ratio)  associated  with  (1.1).  Further,  let  C  and  G 

#"w  /V 

stand  for  the  right  and  left  Cauchy-Green  deformation  tensors,  whence 

C=  FTF  ,  G  =  F  FT  on  ft  .  (1.3' 

Both  £  and  G  are  symmetric,  positive-definite  tensors,  which  have  the 

same  fundamental  scalar  invariants  I.  and  hence  possess  common  positive 
2 

principal  values  x..  ;  consequently,  x^  >0  are  the  principal  stretches 
of  the  deformation  at  hand.  Thus, 


12  =  ^-^tr£-)2  -  tr(£2)]  =  x?x2  +  x|a3  +  x|x^  ,  „ 

13  =  det  C  =  =  X^x^X^  . 


(1.4) 


Since  the  material  is  assumed  to  be  incompressible,  only  locally 
volume-preserving  deformations  are  admissible.  Therefore, 


J  —  X^  X 2 x 3  =  1 ,  =  1  on  ft  . 


(1.5) 


Let  a,  defined  on  ft,  be  the  nominal  (Piola)  stress-tensor  field 

^If  M  is  a  three-by- three  matrix  with  elements  M-h,  we  alternatively 
write''”  [M-|j]  in  place  of  M;  6-jj  is  the  Kronecker-delta. 

2 

A  superscript  T  will  always  indicate  transposition. 
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accompanying  the  deformation.  The  equilibrium  balance  of  linear  and  an¬ 
gular  momentum,  in  the  absence  of  body  forces,  then  demands  that 

div2,=  £>  SL^’lZ  or  °ij,j-°.  °ikFjk  =  Fi*aj*  on  R-  (1‘6)1 

Further,  suppose  £  is  the  true  (Cauchy)  stress-tensor  field  on  ft*.  Ac¬ 
cordingly, 

T(y(x))  =  a(x)FT(x)  V  x  6 ft  .  (1.7) 

Turning  now  to  the  governing  constitutive  relations,  we  call  W  the  elas¬ 
tic  potential  of  the  hyperelastic  material  under  consideration.  The 
scalar-valued  response  function  W,  which  represents  the  strain-energy 
density  per  unit  undeformed  volume,  is  taken  to  be  defined  and  at  least 
twice  continuously  differentiable  on  the  set  £  of  all  nonsingular  second- 
order  tensors.  The  appropriate  constitutive  Taw  may  then  be  written  as 

a- W  -  "CT  or  «,j  *  »W(£)/SF1J  -  pFjj  ,  (1.8)2 3 

in  which  p,  for  the  present  assumed  to  be  continuously  differentiable  on 
ft,  stands  for  the  arbitrary  pressure  field  needed  to  accommodate  the  con- 
straint  of  incompressibility  (1.5). 

The  elastic  potential  W  is  subject  to  the  requirement  of  material 

Subscripts  preceded  by  a  comma  indicate  partial  differentiation  with 
respect  to  the  corresponding  material  Cartesian  coordinate. 

2  -T 

Here  and  in  the  sequel,  F  denotes  the  transposed  inverse  of  the  ten¬ 
sor  F,  while  F“J-(F"T)ij. 

3 

We  emphasize  that  although  this  constraint  restricts  the  argument  of 
W(£)  to  unimodular  tensors,  the  particular  manner  in  which  the  domain  of 
definition  of  W  is  extended  to  the  set  £  of  all  nonsingular  tensors 
affects  merely  the  pressure  p,  and  is  therefore  irrelevant. 
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frame  indifference,  which  demands  that 

W(F)  =  W(£F)  V  (Q,F)  6  OxX  ,  (1.9) 

where  0  is  the  set  of  all  proper  orthogonal  second-order  tensors.  More¬ 
over,  (1.8)  and  (1.9)  are  found  to  imply  the  second  of  (1.6)  and  hence  the 
symnetry  of  the  true  stress-tensor  field  t. 

Substituting  from  (1.8)  into  the  first  of  (1.6),  appealing  to  (1.2), 
and  recalling  that  FT ]  .  =  0  for  a  unimodular  deformation-gradient  field, 
one  is  led  to  the  displacement-equations  of  equilibrium.  Adjoining  to  the 
latter  the  incompressibility  condition  (1.5),  one  arrives  at  the  system  of 
partial  differential  equations 

Sjkt<£K,trp.jFj!-0-  J*det£=1  R.  <'■'<» 

provided  are  the  components  of  the  fourth-order  tensor  defined 

by 

=  W£> *  A(£)/8Fij9Fi«  •  c-11' 

In  the  case  of  material  isotropy,  the  strain-energy  density  W(£) 
involves  F  only  through  the  invariants  I,  and  I,  of  C=FTF.  Thus, 

W(F)  =  W(Ii(F),!2(F))  .  (1.12) 

From  (1.3)  and  (1.4)  follow 


3lo 


ZGuFtJ 


0.13) 


For  future  convenience  we  adopt  the  abridged  notation 
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“a  =  3W(  I-j ,  I2)/3lot  .  Wa6=3ZW(Ii,I2)/3la3l{ 


0.14)' 


Equations  (1.8),  (1.12),  (1.13),  (1.14)  then  yield  the  stress-deformation 
relations  for  isotropic  incompressible  hyperelastic  solids  in  the  form 


-1 


2<  v 'I  VFir  pFji 


or  o  »  2<W,  ♦  I,U2)F  -  2W.& F  -  pF 


-T 


*  (1.15) 

> 


On  the  other  hand,  on  account  of  (1.3),  (1.7),  (1.15),  the  true  stress 
field  £  obeys  the  constitutive  relation 

r.2(ii1M,U2>G-2«.,G2-rtL  ,  (1.16) 

where  1_  stands  for  the  idem  tensor  with  the  components  6^. 

According  to  (1.16),  the  principal  axes  of  ^  and  G  coincide; 
further,  if  designates  the  principal  true  stress  associated  with  the 
principal  stretch  A.. ,  one  evidently  has 

x.j  =  -p  +  2a?[W^  +  (I-j  -  A?)^]  (no  sum)  .  (1.17) 


Later  on,  we  shall  need  to  refer  to  the  Baker-Ericksen  inequalities, 
which  require  that 

(Ti  "  TjX?i  "  Aj)  >  0  lf  Xi^Aj  (no  sum)  ,  (1.18) 

and  thus  postulate  that  the  greater  true  principal  stress  occurs  always 
In  the  direction  of  the  larger  principal  stretch.  Because  of  (1.17)  and 
(1.5),  the  Inequality  (1.18)  is  equivalent  to 

Throughout  this  paper  Greek  subscripts  have  the  range  (1,2). 
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^  +  (x.xj)'2W2>°  if  .  (1.19) 

With  a  view  toward  defining  the  notion  of  ordinary  el  1 ipticity  rele¬ 
vant  to  the  system  of  partial  differential  equations  (l.iO)l  we  consider 
a  surface  W'  lying  wholly  within  ft  that  admits  a  twice  continuously 
differentiable  and  one-to-one  parameterization 

x  =  x(crC2)  *  Ult52)  6  Z  ,  (1.20) 

where  z  is  a  region  of  the  parameter-plane  and  (c1  ,S2)  are  orthogonal 
curvilinear  coordinates  on  Let  P  be  an  arbitrarily  chosen,  fixed 
point  on  r/.  Then,  within  a  three-dimensional  neighborhood  of  P,  we  may 
introduce  orthogonal  curvilinear  coordinates  (c^.S2»c),  such  that 

^=x(c1,C2.0  =  x(c1,C2)  +  cn(c1,52)  ,  (1.21) 

in  which  n^  is  the  unit  normal  vector  of  ^  and  |c|  is  the  perpendic¬ 
ular  distance  from  of  a  point  with  position  vector  x.  The  mapping 
(1.21)  is  locally  one-to-one;  let  its  inverse  be  given  by 

C  «i  (x),  C  =  dx)  ,  (1.22) 

&  Cl 

where  £a  and  £  are  defined  and  twice  continuously  differentiable  in 
a  neighborhood  of  P. 

We  now  weaken  the  original  smoothness  requirements  on  u^  and  p. 

Thus,  we  assume  that  (u,p)  is  a  solution  of  (1.10)  with  u^  continuously 
differentiable  and  merely  piecewise  twice  continuously  differentiable  and 

What  follows  is  at  once  an  adaptation  to  Incompressible  bodies  of  the 
analysis  in  Section  1  of  [1]  and  a  generalization  to  three  dimensions  of 
the  development  in  Section  3.1  of  [6]. 
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p  continuous  but  only  piecewise  continuously  differentiable  on  ft.  We 
shall  call  such  a  (u,.p)  a  "relaxed  solution"  of  (1.10).  On  setting 


u(c1,C2,t)  ,  p(c1  ,e2,c)  =  P(x(51  ,c2,t))  , 

one  confirms  with  the  aid  of  (1.22)  that 


uk,£j~  35a  a,*j  3 c 

2-  2- 
3  \  .  3  Ul 

r  .  4-  - - - F  F  .  +  - — 

,ij  3 C o 3 C  ct,£  B,J  3C3C 
p  ci  a 

9- 

2— 

3  U.  A 

3  U. 
k  a  a 

4-  *  ~ 

3Sa3C  S*S,j 

2  ^  i  » 

dC  ’ 

£  ,• 


(1.23) 


(1.24) 


p  ,-!§-£  i  • 

,J  3^a  3£  *3 


On  account  of  the  assumed  smoothness  of  u,,p  and  of  the  mapping 

(1.21),  the  first  and  second-order  partial  derivatives  of  u^  and  the 

—  2—2 

first-order  partial  derivatives  of  p,  except  possibly  3  uk/3£  and 
3p/3c ,  are  continuous  in  the  neighborhood  of  P  under  consideration; 
further,  the  latter  two  derivatives  may  at  most  have  finite  jump  discon¬ 
tinuities  across  t/.  Moreover,  (1.24)  and  the  smoothness  of  the  inverse 
mapping  (1.22)  give 

Iukaj.]J  =  lE32uk/3c2]]c j  ,  CpjU3  ITap/acIcj  on  J  ,  (1.25) 

where  [[hj  denotes  the  jump  of  a  function  h  across  y. 

The  second  of  (1.10)  implies  J  .  =  0,  which  -  because  of  (1.2),  (1.5), 

9  J 

and  the  last  of  (1.13)  -  is  equivalent  to 


FIkuk,«j=0  8 


(1.26) 
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Consequently,  in  view  of  (1.25),  (1.26)  and  the  continuity  of  F  and  of 
cijk*’  e<luations  (1  -10)  yield 


j  =  0  on  . 


(1.27) 


Also,  noting  that  v£/|v£|  on  >/  coincides  with  the  unit  normal  vector 
n  of  ^  and  defining 


Qik(£;£)  =  cijk*(£>jn*  v  (F,n)€Xx“U  , 
v=[[32u/3c21],  g=  Iver’lap/ac]!  on 


(1.28)1 


where  h.  is  the  set  of  all  unit  vectors,  we  infer  that 

CV  2.’  v-(F"Tn)  =  °  on*/.  (1.29) 

Clearly,  (1.29)  constitute  four  linear  homogeneous  algebraic  equations 
in  the  jumps  (v^  ,q),  which  admit  only  the  trivial  solution  ^  =  0,  q  =  0 
if  and  only  if  the  determinant  of  their  coefficient  matrix  fails  to  vanish. 

The  system  (1.10)  is  said  to  be  elliptic  at  a  relaxed  solution  (u,»p) 
and  at  a  point  x.€ft  if  and  only  if  u,  Is  twice  continuously  differentiable 
and  p  is  continuously  differentiable  at  £.  Accordingly,  when  (1.10)  is 
elliptic  at  (jlj ,p)  and  x,  there  does  not  exist  any  surface  (of  the  req¬ 
uisite  smoothness)  through  the  point  x  across  which  the  "second  normal 

2—2  — 
derivative"  3  u/3?  or  the  "normal  derivative"  3p/3c  is  discontinuous. 

It  is  evident  from  (1.29)  that  (1.10)  is  elliptic  at  a  relaxed  solution 

T - 

Here  Q-ji,  are  the  components  of  the  acoustic  tensor  Q,  which  is  syrnnetric 
because  of  (1 .11 ). 
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(u^.p)  and  at  a  point  £  if  and  only  if 


A(F;n )  =  det 


Qn 

*12 

Q13 

-"l" 

*21 

q22 

Q23 

"n2 

't  0  V  n  6  "U  , 

Q31 

Q32 

*33 

"n3 

."l 

n2 

n3 

0 

(1.30) 


a-ates;.  a=r  n»  z=i+va(>o 


or  equivalently 


A^O)=Ieijkctmn"i"tQjni0kn?<0  v 


Sr&Dn.)*  i=Z'  a-  * 


(1.31) 


where  stands  for  the  components  of  the  three-dimensional  alternator. 

If  (1.31)  is  violated  for  some  unit  vector  n,then  n  is  normal  to  a  ma- 

terial  characteristic  surface  in  R  through  x.  These  characteristic 

2—  2 

surfaces  are  the  only  possible  carriers  of  discontinuities  in  a  u/ac 
or  ap /ac,  and  ordinary  ellipticity  precludes  the  existence  of  real  char¬ 
acteristic  surfaces. 

We  turn  next  to  the  definition  of  strong  ellipticity^  appropriate  to 
the  system  of  partial  differential  equations  (1.10).  For  this  purpose  we 
consider  an  infinitesimal  plane  displacement-pressure  wave  superposed  upon 
a  finite  homogeneous  deformation  of  an  incompressible  :,yperelastic  body 
occupying  the  entire  three-dimensional  space  t.  Thus,  if  t  denotes 


This  concept  is  essentially  identical  with  the  notion  of  Hadamard  sta¬ 
bility.  What  follows  is  parallel  to  the  analysis  in  Section  3  of  [12], 
which  is  however  confined  to  isotropic  materials. 


0  o 

the  time,  while  F  and  p  designate  the  constant  deformation-gradient 
and  pressure  fields  associated  with  the  homogeneous  pre-deformation,  we 
set 

y(x,t)  =  £x  +  w(x,t),  p(x,t)  =  p+ r(x,t)  V  (x,t)  € e x  (-00  ,00)  ,  (1.32) 

regard  |vw|,  as  well  as  |r|,  small  compared  to  unity,  and  take 

w(x,t)  =  a9(t  •  Fx-ct),  r(x,t)  =  *(«,-  F  x  -  ct)  .  (1.33) 


Here  a  and  1  are  constant  unit  vectors,  determining  the  direction  of 
motion  and  the  direction  of  propagation  of  the  plane  wave  at  hand,  whereas 
Cj^O  is  its  speed  of  propagation.  Further,  cp  and  ip  are  scalar-valued 
functions,  the  first  of  which  is  twice  -  the  second  once  -  continuously 
differentiable  on  (-00,00).  Finally,  we  take  for  granted  that  the  deriv¬ 
atives  cp"  and  <|/  fail  to  vanish  identically.  From  (1.32)  and  the  in¬ 
compressibility  requirement,  one  has 


F  =  vy  =  F  +  vw ,  J  *  det  F  =  1  on  ex  (-00,00),  det  F  =  1  .  (1.34) 


The  time-dependent  nominal  stresses  induced  by  the  motion  (1.32), 
(1.33)  follow  from  the  constitutive  relation  (1.8)  and  are  now  subject 
to  the  stress  equations  of  motion 

2  2 

0..  .  =  p3  u,/3t  on  ex  (-00,00)  ,  (1.35) 

IJ  »J  1 

provided  p>0  is  the  constant  mass  density  and  u  is  the  displacement 
field.  In  view  of  the  first  of  (1.32),  we  see  that 

Ujf&t)  “ >i (X.,t)  -  xi  =  (Fu  -  «ij)xj  ♦w1(x,t)  .  (1.36) 


Recalling  the  identities 


f 
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det£4e1jkWF1pF4qFkr>  (<tet&Fj!4e1kt'jpqVtq  V  •  ('’37) 

we  substitute  from  (1.34)  into  (1.37)  and,  upon  a  linearization  with  re¬ 
spect  to  vw,  arrive  at 


J  =  detF-detF  +  l£l.jk£pqJlFipFjqwk>1  =  l+F^wk>Jl  ,  > 


Fj1~Fj1  +  em£jpqFkpVq  • 


y  d.38) 


V 

Further,  expanding  Wp(F)  in  (1.8)  as  a  Taylor  series  in  F  around  F>, 
appealing  to  (1.11),  (1.34),  as  well  as  (1.32),  (1.38),  and  linearizing 
with  respect  to  vw  and  r,  we  obtain 


3W(F) 


iJ  3Fij 


o-PFji+Cijkz^wk,t-Peik£ejpqFkpwt,q-rFji  *  (K39) 

F  *  F 


Now,  (1.35)  and  the  second  of  (1.34),  together  with  (1.39),  (1.36),  and 
the  first  of  (1.38),  lead  to  the  linearized  version  of  the  displacement- 
equations  of  motion  and  the  linearized  incompressibility  condition 

°"  ex  (-00,00)  .  (1.40) 

For  the  plane  wave  characterized  by  (1.33),  equations  (1.40)  give 


[Rik(F;t)ak  -  pc2ai  V(4.-  F  x  -  ct)  =  l^'(l  •  F  x  -  ct) 


(1.41) 


in  which 


R^k^F;^  =  Rki  =  * 


(1.42) 
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From  (1.41)  follows 

ak4k  ^  ^~*  ~2-'  ct)  =  bcf  (i*  Fx-ct),  b  =  ,  (1.43) 

and  hence, 

bi=  a*  £.*i=o-  o-44) 


The  system  (1.44)  constitutes  four  linear  homogeneous  algebraic  equations 
in  (a^.b)  and  has  a  nontrivial  solution  if  and  only  if  the  determinant 
of  its  coefficient  matrix  equals  zero.  It  is  easily  confirmed  that  this 
determinant  may  be  written  as 


rpc2 

R1 2 

R13 

”*1 

R21 

R22-p^ 

R23 

~lZ 

R31 

R32 

R33~pc2 

~h 

*1 

*2 

*3 

0 

?eijkepqr*i  VRjq  ”  PC  6jq^Rkr"pC  5kr^  *  (l-45) 

On  expanding  the  right-hand  side  of  (1.45),  one  sees  that  equations  (1.44) 

o 

have  a  nontrivial  solution  (a^.b)  if  and  only  if  pc  and  the  unit  vec¬ 
tor  £  sati sfy  the  secular  equation 

pc  -  (Rk|(  "  Rij*i*j)pC  +7eijkepqr*i*pRjqRkr  =  °*  ♦  0.46)1 

o 

which  is  a  quadratic  equation  in  pc  . 


Upon  specialization  to  an  isotropic  material  and  to  a  pure  homogeneous 
pre-deformation,  (1.46)  reduces  to  equation  (3.12)  of  Sawyers  [12],  who 
sets  aa(pc2H. 


We  show  next  that  the  discriminant  of  (1.46), 


A  (£’£)  -  (Rkk 


E.’  l(hi)  . 


is  non-negative  for  every  nonsingular  £  and  every  unit  vector  £.  To 
this  end  let 


al  =  R22  "  R33’  a2  =  R33  ‘  R11  ’ 


Y1  =  2R23*  y2  =  2R31 *  Y3  =  2R12  ’ 

o  p  2 

and  note  first  that  if  £^  =  £3  =  0,  one  has  a(£;ji)  =  +  y2  *  0.  On  the 

2  2 

other  hand,  if  £^  +  £3>0,  one  verifies  after  a  straightforward  computation 
making  use  of  =  1 ,  that 

a( F ;i)  =  £2  +  a2t2  -  a3£2  -  y^  *2  "  ^  *2*3  +  Y22'3s'l 

+  2(£2  +  £2)-1 (y2*3  -  a2Z1 )£1£2]2 

+  [y-j£]  -  y242  "  Y3a3  +  2(£.-j  +  £3)  (y2£3  -  a2£-]  )^2£31  ^  0  • 


Consequently,  for  every  £€£  and  every  £_€t<  there  exist  two  -  possibly 

o 

coalescent  -  real  values  of  pc  satisfying  the  secular  equation  (1.46). 

We  proceed  now  to  the  relevant  definition  of  strong  ellipticity. 

The  system  (1.10)  is  said  to  be  strongly  elliptic  at  a  solution  (u_,p) 
and  at  a  point  if  and  only  if  a  body  of  the  same  material  occupying 

6,  having  been  subjected  to  a  homogeneous  deformation  with  the  deformation- 

O  e 

gradient  £=  £+  vu^x)  and  to  a  constant  pressure  p  =  p(x) ,  admits  only  in¬ 
finitesimal  plane  displacement-pressure  waves  of  the  form  (1.33)  with  real 
non-zero  propagation  speeds.  Unless  (1.10)  is  strongly  elliptic  at  (u_,r' 
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and  x,  there  is  an  such  that  (1.46)  is  satisfied  for  two  conjugate 
imaginary  values  of  c  or  for  c  =  0.  In  the  first  instance,  there  exist 
superposed  infinitesimal  motions  of  the  form  (1.33)  with 

w(x,t)  =  aexp[i(*  •  Fx-ct)],  i2=-l  , 

which  grow  unbounded  in  time.  When  c  =  0,  in  turn,  one  can  satisfy  the 
linearized  equations  of  motion  with 

w(x,t)  =  a  t  exp(it  •  £ x)  , 

so  that  in  either  event  the  body  is  dynamically  unstable. 

Evidently,  necessary  and  sufficient  in  order  that  (1.10)  be  strongly 

elliptic  at  a  solution  (u,,p)  and  at  a  point  x  is  that  both  real  roots 

2  0 
pc  of  the  secular  equation  (1.46)  be  positive  for  £=  1+ and  for 

every  unit  vector  i.  According  to  (1.46),  this  is  the  case  if  and  only  if 


(1.47) 


(1.48) 


so  that  the  strong-ell ipticity  conditions  (1.47)  may  be  expressed  in  terms 
of  the  acoustic  tensor  Q  as 


(1.49) 
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Finally,  by  comparing  (1.31)  with  the  first  of  (1.49),  we  confirm  that 
strong  ellipticity  implies  ordinary  ellipticity. 

It  is  essential  to  remark  that  conditions  (1.31),  as  well  as  (1.49), 
are  inadequate  as  a  means  for  testing  the  local  ellipticity  of  a  par¬ 
ticular  equilibrium  solution  appropriate  to  a  specific  material  within 
the  class  under  consideration.  For  the  important  subclass  of  isotropic 
materials,  however,  one  may  deduce  from  (1.31)  and  (1.49),  corresponding 
intrinsic  ellipticity  criteria  in  a  form  suitable  for  this  purpose.  It 
is  this  task  that  constitutes  our  main  goal  and  to  which  we  presently 
turn  our  attention. 


- -  ...  .. 
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2.  Explicit  necessary  conditions  for  ordinary  ell ipti city  in  the  case 
of  isotropy. 

Our  current  objective  is  to  deduce  from  (1.31),  which  is  necessary 
and  sufficient  for  ordinary  ellipticity,  a  set  of  explicit  necessary 
ellipticity  conditions  for  the  special  case  of  an  isotropic  material. 
This  set  of  conditions  involves  exclusively  the  local  principal 
stretches,  which  enter  directly,  as  well  as  through  the  first  and  second 
partial  derivatives  of  the  strain-energy  density  with  respect  to  the 
deformation  invariants  or  the  principal  stretches. 

As  a  first  step,  we  derive  an  ellipticity  condition  equivalent  to 
(1.31),  but  analytically  more  amenable.^  Since  [C.  .]=FTF  is  a  sym- 

i  j  >•  <«■ 

2 

metric  matrix  with  the  principal  values  Xi  ,  there  is  an  orthogonal 
matrix  M=  M(F)  ,  such  that 

MTCM  *  D  ,  D. .  =  x?  6,,  (no  sum)  .  (2.1) 

Next,  for  every  nonsingular  F^  and  every  unit  vector  m  ,  let 
N*N(F)  and  H=H(F;m)  be  the  auxiliary  matrices  defined  by 

******  r* 

N(F)  =  mV1  ,  H(F;m)  *  NQ(F;Mm)  NT  ,  (2.2) 

A*  IV  ^  ***  rs* 

in  which  Q  is  the  component  matrix  of  the  acoustic  tensor  introduced 
in  (1.28).  Evidently  H  is  symmetric.  On  setting 

*0 

51  =  MT£,  n=F"Tn,  m=  Nn  ,  H  =  H(F;m)  ,  Q=Q(F;n)  V  (F,n)  €  £*  V 

(2.3) 

What  follows  is  suggested  by  the  development  in  Section  3  of  [1]. 


one  readily  confirms  the  block-matrix  equation 


Appealing  to  (1.30)  and  (2.1)  through  (2.4),  one  sees  after  elementary 
manipulations  that  for  every  unimodular  F  and  every  unit  vector  m. 


A(F;Mm)  =  det 

to* 


'll 

H12 

H1 3 

-s 

'21 

H22 

H23 

-m^ 

,  m  *  D-1m  , 

'31 

H32 

H33 

-m3 

/w  /w 

'l 

m3 

0 

where  M  is  the  orthogonal  and  D  the  diagonal  matrix  appearing  in 
(2.1).  Consequently,  a  necessary  and  sufficient  condition  for  ordinary 
ellipticity,  equivalent  to  (1.31),  is  given  by 


fl(L;S  =  I  EijkGpqr^i^pHjqHkr  *  0  '  •  T  E  V 


Next,  from  (1.3)  and  (1.11)  through  (1.14),  one  finds  that  for  an 


isotropic  material 
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c1Jki^ 


=  2(ft]  +  I1W2)'Sik6jJl  +  4(W2  +  W11  +  21^^+  ^  W22)  Fi jF|cfi. 

•  2Q2<,:ie|:kj*G1k5jltS1kCjt)-4("l2+Ili22><FioGkn,Fin1  f<2'6) 


+  GimFmjFk£^ +  4^22GimFmjGknFn£ 


Further,  (2.2),  (1.28),  (2.6),  (1.3),  (2.1),  and  the  first  of  (2.3)  even¬ 
tually  yield 


H  (F;m)  ■  2CV<Ir  “ktVk^lj  *  4(?  “2  *  “ll  *  2¥l2  +  ‘l"22>Vj 


-  2Vij-4(“l2+Il"22>(m1DjkmktDit"t,"j)t4“22DikmkDjA  •  J 

(2.7) 


Let  ,wi  and  B^-.B  be  defined  by 


8,  ■  “1  *  a;«2  • 


wi  a  Ai  "  Ai  ^8i  +  2Ai^l  "  2I2  “  Ai  '  2Ai2^^ll  +  2Ai  W12  +  Ai  W22^  M2-8^ 


.  2  ,  *  ,  ,  4  , , 


(no  sum)  , 


BU  =  Al2fi203  ’  B22=A22s38l  ’  B33  =  A320102  *  B12  =  02W3  +  A1  0103  ’ 

B21  =  81w3  + A220203  »  B23  =  03W1  +  A22s201  ’  B32  =  02W1  +  A3  0301  ’ 

B31  =  01w2  +  A320302  *  B13=03W2  +  A1  0102  ’ 

o  ?  _?  2  -22  -42  -42  -42  „  2 

B  =  4X^8^ +  4X2  82  +  4x3  83  "  xi  *V  A2  w2  -  *3  w3  +  2a3w1w2 

2  2 
+  2x-j  w2w^  +  2^2  w^w-i 


(2.9) 


We  from  (2.9),  (2.8),  (1.14),  (1.4)  that  B.  .  and  B  are  solely 
functions  of  the  A^. 
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Upon  a  lengthy  computation  involving  (2.5),  (2.7)  through  (2.9),  (2.1), 
and  (1.5),  one  verifies  that  for  every  unimodular  F  and  every  unit  vec- 
tor  m  , 

A(F;Mm)  =  E(x_;z)  i  4Bijz?zj  +  Bz^^  ,  zk  =  m*  .  (2.10)1 

Suppose  now  that  (1.10)  is  elliptic  at  (u,p)  and  x  .  If 
x^  =  Xi(x)>0  ,  with  x-jX2X3=  1*  are  the  corresP°nciln9  local  principal 
stretches,  one  thus  has 


E(X;z)  i  0  V  Z  €  A  ,  (2.11) 

A#  fV  *** 

where 

A  =  { z  |  i  0  ,  z1  +  z2  +  z3=l}  .  (2.12) 

Evidently,  A  is  the  bounded  and  closed  plane  region  in  £-space  whose 
boundary  is  the  equilateral  triangle  (Figure  1)  with  vertices  at 

5,  =  (1,0,0)  ,  ^  =  (0,1,0)  ,  £3  =  (0,0,1)  .  (2.13) 

In  view  of  (2.11)  and  the  continuity  of  E(*;z)  in  z  ,  this  function 

/w 

is  of  one  sign  on  A  ,  so  that  in  particular. 


E(x;CJ  E(X;C.)  >  0  ,  E(X;c„)  E(X;cJ  >  0  . 

r**  rv  I  rw  /m  a*  /vJ 

Consequently,  by  virtue  of  (2.10),  (2.9), 

3^83  >0  ,  ^2^3  >  ^  *  (2.14) 


and 


^Here  and  in  the  sequel,  £  and  z  stand  for  the  triplets  of  real  num¬ 
bers  (Xj,x2,x3)  and  (z~,z2,z3)7  respectively. 


1 
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E( A;z) >  0  V  z  € A  .  (2.15) 

In  order  to  deduce  additional  necessary  conditions  for  ordinary 
ellipticity,  we  require  the  following 

LEMMA.  Let  a,b,c,d  be  real  constants,  with  a>0  and  d>0. 

Then 

as3  + bs2t  + cst2  +  dt3  >  0  for  all  s>_0,  t>_0,s  +  t=l,  (2.16) 
if  and  only  if 

either  4ac3  +  4db3  -  b2c2  -  18abcd  +  27a2d2  >  0  or  b>0,c>0.  (2.17) 

To  prove  this  lerrma,  note  first  that  (2.16)  implies 

0(c)  =  ae3  +  be2  +  cc  +  d  >  0  for  all  C  >.  0  . 

Necessary  and  sufficient  in  order  that  the  cubic  polynomial  e(c)  have 
three  real  zeros  is  that 

4ac3  +  4db3- b2c2- 18abcd  +  27a2d2  <  0  .  (2.18)1 

Accordingly,  (2.18)  together  with  (2.16)  imply  the  existence  of  real 

numbers  £.  (i  s  1,2,3)  and  T  (a*  1,2)  ,  such  that 
1  a 

9(5i)  =  0  ,  9'(ca)*0,  5,  <5]i52<C2S53<0  . 


Thus 


e'U)  *  3ar  +  2bc  +  c  =  3a(c  -  C-,)(C 


c2) 


^See,  for  example,  Dickson  [14],  p.  47.  The  truth  of  this  assertion 
requires  merely  that  a,b,c,d  be  real  with  a^O  . 


m iigMssiausfe 
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and  therefore  (2.16)  and  the  negation  of  the  first  inequality  in  (2.17) 
necessitate 


b  =  -|a(c1  +  C2) 


c  =  3a  5^2  >  0  • 


Hence  (2.16)  implies  (2.17).  Conversely,  b>0  and  c>0  at  once  ass  .ire 
that  (2.16)  holds.  On  the  other  hand,  (2.16)  also  follows  from  the  first 
of  (2.17).  Indeed,  the  latter  mandates  the  existence  of  but  one  real 
zero  of  0(C);  this  zero  is  negative  because  a>0  and  d>0,  so  that 
0(C)  >  0  for  all  CiO,  which  is  easily  seen  to  give  (2.16).  This  com¬ 
pletes  the  proof. 

The  inequality  (2.15)  in  particular  requires  E(*;z)>0  on  the 
side  of  the  triangular  boundary  of  A  with  z^  =  0,  whence  from  (2.10), 

B22Z2  +  823Z2Z3  +  B322223  +  B33Z3>  0  for  a11  z2~°’  z3~0,  z2  +  Z3  =  1  ' 

(2.19) 


The  foregoing  lemma  thus  entitles  us  to  conclude  that 


either 

or 


4B22B32  +  4B33B23  "  B23B32  "  18B22B23B32B33  +  27B22B33  >  0  I 


B22  >  0  ,  >  0  • 


(2.20) 


Bearing  in  mind  (2.14),  we  now  set 

pi  =  S^/X - e-| B263  >  0  ,  u.  =  w./x.jg.  (no  sums)  , 


(2.21) 


and,  with  the  aid  of  (2.9),  (2.21),  infer  the  identity 
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4B22B32  +  4B33B23  '  B23B32  '  18B22B23B32B33  +  27B22B33 


-  .4.4*2  2  r  t  — 1  ,  “1  \*i2/2  *\ 

—  A.-J  S-j  $2  63  Lw-j  -  (p2p3  +  p2  p3 “  ^ ) 


(2.22) 


Further,  on  account  of  (2.22),  (2.9),  (2.21),  we  draw  from  (2.20)  that 


either  u>^  <  4 


or  u>.j  +  P2P3I  >  0  ,  u-j  +  p3  >  0 


(2.23) 


Since  the  second  alternative  in  (2.23)  implies 
u>i  >  -min(p2p3^ ,p2^ P3) i  -1  * 

one  has  u-j  +  2  >  0  .  Thus  and  from  parallel  arguments  applied  to  the 
edges  of  A  lying  in  the  planes  z2  *  0  and  z3=0,  one  arrives  at 

<ui  +  2  >  0  ,  (2.24) 


as  necessary  conditions  of  ordinary  ell ipti city. 

With  a  view  towards  applying  (2.15)  to  certain  points  in  the 
i nterl or  a  of  the  triangular  region  A  ,  we  now  suppose  that  z^  >  0 
and  define  functions  vi  through 


v2 

v3 


P1  (“l  +  X3X21z2Z31  +  X2X31z3Z2^  +  °2X3X1 1 Z1  Z3] 
P2^2+X1X3  Z3Z1  +X3X1  zlz3]  ^  +  °3X1X2  Z2Z1] 
p3((i)3+x2xi  z^2  Z2Z1 1^  +  P1X2X31  Z3Z21 


+  p  3X2  X1 
+  0 1 X  3X2 
+  P2X1X3 


zlz2  ’ 

z2Z31,V(2.25) 


After  some  tedious  algebra  based  on  (2.10),  (2.9)  and  involving  (2.25), 
(2.21 ),  one  obtains 
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E(*.;z)  =  *  e3zlz2z3  (v?  +  v2  +  v3'  2v1v2"  2v2v3‘  2v3v1  > 

for  >  0  .  (2.26) 

Consequently,  (2.15)  yields 

V1  +  v2  +  v3  "  2v1v2  "  2v2v3  “  2v3v1  <  0  on  ^  •  (2.27) 

and  thus  from  (2.25),  (2.24),  (2.21)  one  infers  that  (2.15)  necessitates 
v3  <  V1  +  v2 +  2^?2  *  v.  >  0  on  A  .  (2.28) 


For  our  present  purpose  it  is  expedient  to  introduce  the  transforma¬ 
tion 

n  =  n  (z;A)  =  X-A_1z  zZ1  V  z  6  A  (no  sum)  ,  (2.29) 

(X  Qt  m  J  Q  Cl  J 

which  is  a  mapping,  depending  parametrically  on  x^  ,  of  A  into  the  open 
first  quadrant 

n  *  { (r»i  >r»2 )  I  n-j  »  0  ,  ng  >  0  } 

of  the  (n-j ,n2)~plane.  Inverting  (2.29)  -  keeping  in  mind  that 

O 

z1+z2  +  z3=l  on  A  -one  is  led  to 

_1 

za  =  Za^nl  *n2*^  ~  ^aria^lr'l  +  ^2^2  +  *3^  (no  Sum)  , 

(2.30) 

—  z3(n-| » ng  +  ^2^2  ^  3  ^  ^  ( n  -j  »n2 )  6  n  . 


Thus,  (2.29)  in  fact  constitutes  a  one-to-one  mapping  of  A  onto  n  . 
Next,  let  functions  f  and  g  be  defined  by 
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"f(ni  *n2lw-j  »p-j  »P2»p3)  =  P]  (“i  +  ^2  +  ^2^  +  °2nl  +  p 3nl n2 ^  * 


9 ( n-j  *02 »^2 * P "]  ,p2,p3^  ™  plwl  ^  P2W2  ^  *"p2^  1  ^p1^2 


>(2.31) 


1  /9 

+  2[f  ( n-j  ,n2  *i»>1  »p  1  ,p£ 02 »n-j » <^2 *?2 »p  1  »p^ )]  V  (01*02)  ^ n  J 


in  which  w  ,p.  are  given  by  (2.21),  (2.8)  and  ultimately1  depend  solely 
on  £  .  It  is  clear  from  (2.31),  (2.29),  (2.25),  and  the  second  of 
(2.28)  that 


f(n1(z,x)  ,  02 ( z ; x ) \ (*>1 , p  1 , P2 * p 3 )  —  Vi  >0  , 
f  ( n2 (^-’i.)  *  ^1  *w2  ,p2 *p  1 ,p 3 ^  =  v2  >  ®  V  z  £ 


A  . 


(2.32) 


Further,  as  a  consequence  of  (2.25)  and  (2.30)  through  (2.32),  the  first 
of  (2.28)  is  equivalent  to 


p^**^  <  9(01  tOg*^  »“>2»p']  »p2,p3^  ^  (01*02)  6  n 


(2.33) 


At  this  stage  we  derive  from  (2.33)  an  additional  restriction  on 
»p 1  by  choos 1  ng  01  and  02  so  as  to  mi m rm ze  g(oi *02*^1 *^2 *pi *p2,p3^ 
on  n.  Thus,  holding  «  ,p^  fixed,  we  take  (01  ,o2)  =  (iT,  ,o2)  .  where 


3g/3o„ 


=  0  . 


(2.34) 


( Oi *02 ) 


Substituting  from  (2.31)  into  (2.34),  we  are  led  to  the  unique  solution 


T 


Recall  from  (1.4)  that  Ia  ,  and  hence  W  (Ii ,I2),  tyjgOi.M,  are  ex¬ 
pressible  in  terms  of  the  principal  stretches. 
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n,  =  <0  +k)_1  ,  n2  = 
and  (2.31)  now  gives 
f(nr?r2  »(Di  »P2 

g(hj  ,tt2  $o>2 » P i 


(l+<)  ^  ,  k  =  [p^  (u-|  +  2)/q^2  +  2)]^^ 


2 

»p^)  ”  k  f  ( ri 2 » *^2 >P 2 *P  1  »p 3 ) » 


_  _  ^ 

,p2,p3^  =  ^p3  +  C^p-j  (**»i  +  2)  +  /p2(w2  +  2)]* 


,  ( 2 . 35 ) 1 


|  (2.36) 


Finally,  from  (2.33),  in  view  of  the  second  of  (2.36),  follows 


03(0)3-2)  <  [/p.|  (o)i  +  2)  +  /p2(u>2  +  2)]2  .  (2.37) 

Proceeding  similarly  from  the  two  cyclic  permutations  of  (2.28),  one 
arrives  at 

p-j  (u>-j  -  2 )  <  [•Gjtzpnj  +  /pjtwj + ?)]^  ,  1 

f  (2-38) 

p2^‘°2”2^  <  ^p3^“3 +  +  ^p-j  (o)-j  +  2)]2  .  J 


We  now  summarize  the  results  established  in  this  section.  To  begin 
with,  equations  (2.10),  (2.15),  and  the  first  of  (2.3)  justify  the  asser¬ 
tion:  for  the  special  case  of  isotropy,  the  system  of  partial  differen¬ 
tial  equations  (1.10)  is.  elliptic  in  the  ordinary  sense  at  a  relaxed 
solution  (u,p)  and  at  a  point  x  only  if 


Equation  (2.34)  merely  characterizes  (tT],"^)  as  a  stationary  point 
of  g(n-| *“2»pl »p2»p3^’  that  this  choice  of  (m»n2)  actually 
minimizes  g  on  n  for  fixed  A^-as  is  shown  in,  and  essential  to,  the 
analysis  in  Section  3  —  is  irrelevant  at  present.  Indeed,  as  far  as 
our  immediate  objective  is  concerned,  but  for  lack  of  motivation, 
(tTi.tTo)  could  have  been  defined  directly  by  (2.35),  without  recourse 
to  (2.34). 
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A(F;n)  >  0  y  n^,'  F  =  1  +  7u(x)  , 

^  rw  #v  ^  /v 


(2.39) 


where  A(F;n)  is  the  determinant  defined  by  (1.30)  and  involves  the 
acoustic  tensor  introduced  in  (1.28).  Moreover,  in  view  of  (2.14), 
(2.24),  (2.37),  (2.38),  a  set  of  explicit  necessary  conditions  for 


in  the  circumstances  under  consideration  is  supplied 


by  the  inequalities 


Here  p^.w-.e^  are  accounted  for  in  (2.21),  (2.8),  (1.14)  and  are  given 
by 

p.,-  =  8- j/*i6i02s3  ’  “i  s  Wf/Ai8i  (no  sums)  •  (2.41) 

8i  =  Wi(Ii>I2)  +  xl  W2(ITI2)  ’ 
wi  *  ^?(In  -  )6i  +  2xf(I^  -  2I2  -  x^-  2xT2)[W11(I1  ,I2) 

+  2 X i  U'j 2 ( ^ ^ 2 ^  Ai  ^22^1**2^  sum)  , 


where  x^»x^(x)  are  the  local  principal  stretches  associated  with  the 
displacement  field  u^,  while  ^(Ij,^)  and  W  (Ij.Ig)  are 


in  which 
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( 1 7  (^|  * ^2^  *  ^2^1  **2^  * 


w  a  -  w  fl (1 1 ( x, , ) 1 1« ( x, i x« ) )  . 


>  (2.47) 


aB  oiB'  1t'a]*a2/,2v1*2 
By  means  of  (2.42),  (2.43),  (2.46),  and  (2.47),  we  can  verify  that 

,  l/,2_,2,-l  * 


B-j  =  B-j  (X-j  ,^2)  =  2  (X2-X3)  Xgl^Xj  tXg) 


-  0  9  1  1  n  „ 

w-j  -  W-|  (x^  jXg)  -  -(X2  -  X3)  X2  ^2  (X1  »X2)  "*’2’  *3  ^22  ( *^2)  »  *3  =  (XiXg) 


1  ,-2* 


(X2  ^  X3)  , 


82  s  S2( X1  *^2)  ”  2  (^i  “  X3)  Xi W]  ( x^  ^2)  , 


w2  ~  w£(  Xi  »X2)  -  -(X-j  -  x^)  X-J  W-j(x.|,X2)  +  2  X3  Wii(x-j,X2)  ,  x^=  (x^) 

( X3  f  X-j  )  , 

83  -  83(x^  ^2)  =  2"  (X-J  -  X2)  [Xi (Xj  .Xg)  -  XgW2 ( X^ , X2 )3  » 

*#3 -  w^( x^  jXg)  -  (Xi  -  Xg)  C -X 2  ^2 ( ^2)  ■  X]  W-j  (x-j  1X2)3 

+  2^2  W-j j  (Xj , X2 )  +  Xj  W22 ( X-j  > X2 )  -  2x<j  X2  W-j 2 ( X-j » X2 ) ]  (Xj  j*  X2) .  J 


(2.48) 

Equations  (2.48)  exclude  the  case  of  coalescent  principal  stretches. 
In  order  to  derive  representations  for  6j,w^  applicable  when  two  of 
the  Xj  are  equal,  one  makes  use  of  (2.44), (2.43)  to  infer 

W ( x^ , X2 )  3  W(X2»Xj )  ~  W^.Xj)  =  W^.x-j )  ,  Xj  ®  ( x-j X2 )  , 


whence 


( ^1  =  ^2 (^2*^1  ^  *  W11  =  ^22 ^2*^1  ^  ’  ^12^1  ’^2^  =  ^1 2  ^2*  ^1  ^ »(®) 

^2^2 ( ^1  **2)  =  ^2*^1  ^2*^3^  "  ^3^2 ( ^2 * ^3 ^  *  > 

^1^1  ^1  ’*2^  =  “*3^1  ^3»^] )  ^l^2^3*^1 )  * 


2  * 

*2^22^1  ’*2^  — 

2X3W2 ( x2 •  ^3)  "*"  *2  ^11  ^2*^3^  *  *3  ^22^2’*3^  ”  ^A-j  ^12^2* ^3^ 
A-|  W]  1  ( x -j  ,*2 )  = 

2A3W1  (^3»^i )  +  ^3  ^n  ( A3»X^ )  +  A-j  ^22(^3*  A-| )  -  2A2  ^1 2(^3*^! ) 


*3  =  ( A-j  ^2 ) 


(2.49) 


Substitution  from  (b)  of  (2.49)  into  (2.48)  gives 


"it  "It  "jc  *Ar 

6-j  (A-j  ,A2)  -  8j(  A2»A2)  »  ^"j  (A-j  ’Ag)  _  w3(*2’*3^  ’ 


^2 (^1  * A2 )  =  ® 3 ( A 3 »A]  )  ,  W2 ( A'j  » A2 )  =  ( A^  1 A-j  )  ,  A^  =  (^1  Ag) 


(2.50) 


Further,  on  account  of  the  first  two  of  (a)  in  (2.49),  an  obvious  limit 
★  ★ 

process  applied  to  63, w3  in  (2.48)  yields 
63(A,A)  =  1  A2w3(A,A)  =|  [a’1H1(A,A)+H11(A,A)  -  W^U.A)]  (a  >  0).  (2.51) 

Consequently,  if  one  sets 


S  ( A 1  ,  A  2  )= 


\  [a”1^  (a, a)  +  W-j  1  (a, a)  -  ^1 2  (  a  ,  a  )  ]  for  A1  =  a2=  a  , 
J  ( A 1  ~  A 2 )  [A-|W-]  (A^  ^2)  -  A2W2 (A^.Ag)]  for  A 3  ^  A 2  , 


w ( X i ,X2)  = 


I  x‘2[x"1W1(x,x)  +  W11(x,a)  -l*12(x,x)]  for  X1  =  x2  =  X, 


s^( X 1  —  X 2 )  C^2  W2(Xi  ,X2)  -  X Wi  ( X-j  ,X2 )] 

+  ?  ^X2  Wll(xrX2)  +  Xl  ^22^X1  ,X2^  ~  2xl  X2  ^12^X1  ,X2^ 

for  A]  t  x2> 


(2.52) 


then  (2.48),  (2.50),  (2.51)  justify 


S^j  b(x2»A3)  >  62“6(^2’^]^  *  e3  “  e(x-|  »X2 ) 


w1  w(x2,x3)  *  w2  -  w(x3,  A-j )  ,  w3  -  w(x^ ,x2)  , 


(2.53) 


which  accomnodate  all  x^  >  0  .  Finally,  from  the  second  and  fourth  func¬ 
tional  relations  in  (b)  of  (2.49)  and  the  second  in  (a),  one  finds  that 
B.w  may  more  conveniently  be  written  as 


1  *  _2 

'’g'  «1 1  (a,X  )  for  A]  =  X2=  A  , 


0Ur*2W 

\  1  p  p  ★  1 

x2^xl  ”  x2^  ^l  (*2’x3^  for  X1 1  x2,  ^3=(a,jA2)  , 


w(x1 , X2 ) = 


1  x'^^x.x"2)  for  A1  =  x2=  x  , 
•A^U2- A2)'1  W1(x2,x3)+1  a"2W11(X2,X3) 


for  X'j  ^  x2,  X3  =  (x.jA2)  y 


(2.54) 


Equations  (2.53),  (2.54)  comprise  the  desired  intrinsic  representa¬ 
tions  of  B.j  and  w^  in  terms  of  the  xi  and  the  requisite  partial 
* 

derivatives  of  W.  Moreover,  and  w.  now  exhibit  the  cyclic 
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symmetry  inherent  in  (2.53).  Each  of  (a),  (b),  and  (c)  in  the  el  1 iptic- 
ity  conditions  (2.40)  is  carried  into  itself  by  a  cyclic  permutation  of 
the  indices.^  For  the  purpose  of  writing  (2.40)  as  a  set  of  inequalities 
symmetric  with  respect  to  *  we  define  functions  p.w.cf. 

through 


p U-|  *  ^1^2^  ^^2* ^3^  b(X2»x^)  , 

U)(x1  ^2)  -  X^A  2«  (x^  ,X2)/s(x^  .Xj)  i  A3  -  (  A^) 


>  (2.55) 


cp  *j  (x^  -  b(A2»A2)  b(A2,x^  )  j  ^2  * 

X-J  fXo)  s  p(x-|  ,X2)[2  -  U)(x^  .Xg)]  ^  [  *^P  ( X2 » x  2  )cjj  2  (x  2  >x  2 1 


>  (2.56) 


+  ^pTx^TipTJIx^Txp]  »  X3  -  ( x -j x 2 ) 


j 


with  e.w  given  by  (2.54).  Because  of  (2.55),  (2.53),  (2.41),  we  have 


p(x-| =  P2  »  p(X2»X2)  s  p^  ,  p(X2»X-j)  -  P2 

to(x1  ,X2)  —  ^2  *  w(X2fX2)  ~~  <^1  *  u(x2,x^ )  —  ui2  , 

and  thus,  on  appealing  to  (2.53),  (2.56),  conclude  that 

^p-j  (^1  ,^2^  =  ®l®2  *  ^*i  ( ^2 * A3 ^  ~  ®2®3  *  cp]  ( ^2 * X T )  ~  ^3^1  * 

^2 (X] , X2 )  =  <^2  +  2  ,  ^2^2’*3^  ~  *  ^2 ^3,a1  ^  “  u2  +  ^  * 

93(^1  .Xg)  =  p3(2  -  0)3)  +  [ /p  1  (a»^  +2)  +  /P2U2  +  2)  , 


r- 


57) 


1  After  one  cyclic  step,  (a)  of  (2.40)  becomes  8281  >  11  -  >0,  which 

is  equivalent  to  (a). 
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3.  Sufficiency  of  the  conditions  necessary  for  ordinary  eHipticity. 
Explicit  necessary  and  sufficient  conditions  for  strong  elliptic- 
ity  in  the  case  of  isotropy. 

In  this  section  we  show  first  that  the  set  of  necessary  ell ipti city 
conditions  (2.40)  is  also  sufficient  in  order  that  the  system  of  partial 
differential  equations  (1.10)  be  elliptic  in  the  ordinary  sense,  pro¬ 
vided  the  solid  at  hand  is  isotropic.  Thereafter,  we  prove  that  (2.40) 
with  (a)  replaced  by  8i  >  0  are  necessary  and  sufficient  conditions  for 
for  the  strong  ellipticity  of  (1.10).  Finally,  at  the  end  of  the  sec¬ 
tion,  we  note  circumstances  in  which  the  conditions  of  ordinary  and 
strong  ellipticity  coalesce. 

Aiming  at  the  first  of  the  foregoing  three  objectives,  we  now  sup¬ 
pose  (2.40)  holds  and  demonstrate  that  then 

E(x;z)  >0  »ZU  ,  (3.1) 

where  E(a;z)  and  A  are  given  by  (2.10),  (2.9),  (2.8),  (2.12).  It  is 
evident  from  (2.5)  and  (2.10)  that  (3.1),  in  turn,  will  suffice  to 
assure  the  local  ordinary  ellipticity  of  (1.10)  at  the  relaxed  solution 
under  consideration. 

From  (a)  of  (2.40),  in  conjunction  with  (2.9)  and  (2.41),  directly 
follows 

B^.  >  0  (no  sum)  ,  pi  >  0  .  (3.2) 

Further,  since  w-j  +  2  >  0  according  to  (b)  of  (2.40), 

either  -2  <  <  2  or  >.  2  .  (3.3) 


On  the  other  hand,  (3.3)  and  the  second  of  (3.2)  permit  the  inference 
either  -2<  <  2  or  “l  +  P2P3^  >  ^  ’  W1  +  P2^°3  >  ^  ’ 

so  that,  because  of  (a)  in  (2.40), 

either  -x4e|  *  p^1^)]2^  -  4)  >  0  1 

[  (3.4) 

or  ^  1  ®  1  ^ 3 +  P2P3  ^  >  ^  *  *1  S 1 6 2 (^-j  +  P2  p  3 )  >  ®  •  J 

Next,  recalling  (2.22),  (2.41),  and  (2.9),  we  see  that  (3.4)  is  equiva¬ 
lent  to  (2.20).  But  (2.20),  as  a  consequence  of  the  first  of  (3.2)  and 
the  lemma  established  in  Section  2,  implies  (2.19),  which  then  -  by 
virtue  of  (2.10)  —  gives 

E(x;z)  >  0  for  every  z€ A  with  z,  =  0  . 

Proceeding  similarly  from  the  remaining  two  inequalities  in  (b)  of  (2.40), 
one  thus  obtains 

E(x;z)  >  0  Vz€3A  =  A-  A,  (3.5) 

whence  (3.1)  is  true  on  the  edges  of  A  . 

O 

We  have  yet  to  confirm  the  inequality  (3.1)  on  the  i nterior  A  .  As 
a  first  step  in  this  direction,  we  show  presently  that 

g(n-|  »n2»w-|  «^2,01  ,p2’p3^  —  3^1  ^’“l  »“2*pl,p2’p3^  ^  (n^*^)  6A  *  (3.6) 

where  g  is  the  function  introduced  in  (2.31),  while  Tf  is  given  by 
(2.35);  as  before,  n  denotes  the  open  first  quadrant  of  the  (n^.r^)- 
plane.  Accordingly,  we  are  to  prove  that  the  choice  (n-i*^)3  (n"i  .n^)  , 
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which  is  already  known  from  (2.34)  to  render  g  stationary  for  fixed  x  , 
in  fact  minimizes  this  function  on  n  . 

In  order  to  verify  the  above  claim,  we  note  that  the  first  of  (2.31) 
allows  us  to  write 


2-1  2  -1  * 
f  (n-j  .02^1  >Pi  *02*03)  =  ^*1  +  °ln2  ^  ”  n2^  +  °2ni  "*"^3r*lr*2 

2-12  -1 
^(n2»ni  iu)2»P2»0i  »p3)  -  ^2  +  p2nl  +  0^2  +  P3n-|  ^2  * } 


M3. 7) 


provided 


*  =  /p  U  +2)  (no  sum)  .  (3.8) 

a  a  a 

A  lengthy,  but  straightforward,  computation  then  confirms  that 


f(nvn2  ioj-|  ,Pi  ^iPg)  f(n2>n-j  ;w2*<52 *p  1  ’  p3  ^ 


2 

=  Cp3  + 'P-|'<'2  + p-j(l  _  ^  +  p20  -  n-j)] 

+  (p3p i  ni  ^  +  p2p3n2^  +  P1  p2nl 1  n21  ^nl  +  n2  "  1  ^ 

+  °3n/n21  ^2^1  "  nl’i;2^2  +  p]n2^1  ”  n21  ^  ”  n2^2^ 


+  p2n]^2  "  rV  ^  ”  n-j  )^i 


(3.9) 


Since  p .  >  0  and  ^  >  0  by  (a),(b)  of  (2.40),  all  terms  in  the  right- 
1  a 

hand  member  of  (3.9)  are  non-negative,  and 


f(n-j  »H2»<*>i  »p,j  >p2,p3  )  f(n2  » n-j  ;w2»p2»0-|  ) 

2 

i  Cp3  +  i'i’1'2  +  pi  ^  ”  n2^  +  p2^  ’  nl  ^  V  ( ni » n2 )  €  IT 
Combining  this  lower  bound  with  (3.8)  and  the  second  of  (2.31),  we  are 
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1  ed  to 

2 

g(n-|  »n2i“i  »^2»pi  «p2’p3^  —  ^p3  +  (“i  +  2)  +  ^2^2  +  2)] 

V  (nrn2)  €n  .  (3.10) 

Finally,  from  the  second  of  (2.36),  it  is  evident  that  (3.10)  is  equiva¬ 
lent  to  (3.6). 

Next,  from  (3.10)  and  the  last  inequality  in  (c)  of  (2.40)  at  once 
follows  (2.33).  Furthermore,  (a)  and  (b)  of  (2.40)  assure  that  v.>0 

O 

on  a,  if  the  vi  are  the  three  functions  defined  in  (2.25).  But 
(2.33),  with  the  aid  of  (2.31),  (2.29),  (2.32),  and  (2.25),  is  then 
readily  found  to  imply  (2.28).  Parallel  considerations  applied  to  the 
first  two  inequalities  in  (c)  of  (2.40)  result  in 

V1  <  v2  + v3  +  2,/^3  ’  v2<  v3  +  vl +2,/VT  on  (3.11) 

Therefore,  as  v^  is  positive, 

+  >/vJ  ,  ^  +  /vj  on  A  , 

and  thus 

(/\Tj"-  ^vj)2  <  v3  on  A  .  (3.12) 

The  inequalities  (2.28)  and  (3.12),  in  turn,  enable  us  to  conclude  that 

-2/v-|  V2  <  v3  "  V-|  -  v2  <  2/v1  v2  on  A  , 

whence  we  arrive  at  (2.27).  The  latter,  because  of  (2.26)  and  (a)  in 
(2.40),  yields 


E(A;z)  >0  V  zfcX  , 


(3.13) 
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which  -  together  with  (3.5)  -  justifies  our  original  claim  (3.1). 

We  are  now  in  a  position  to  assert  that  the  set  of  inequalities 
(2.40),  or  alternatively  its  symmetric  equivalent  (2.58),  are  both 
necessary  and  sufficient  for  the  ordinary  ell ipti city  of  the  sys tem  of 
partial  differential  equations  (1.10)  at  a  relaxed  solution  (jj , p )  and 
a  point  x  ,  provided  (*)  are  the  local  principal  stretches. 

Our  next  task  is  to  derive  from  (1.49)  an  explicit  set  of  condi¬ 
tions  necessary  and  sufficient  for  the  strong  ellipticity  of  (1.10)  in 
the  special  case  of  isotropy.  To  this  end  we  recall  (2.2),  (2.3)  and 
note  that 

n  =  NTm  ,  Q  =  N'1  H  N"T  ; 

«v  ^  #v  /w  ****** 

hence 

n-n  trQ-n-Qn  *  (NTm)  •  (NTm)  tr(N-^  H  N~T)  -  (NTm)  •  (N  ^Hm)  .  (3.14) 

Further,  from  (2.2),  (2.1),  and  (1.3)  follow 

MT  3  2  1  •  tKjf1  H  N’T)  •  tr(OH)  . 


so  that  (3.14)  gives 


-1. 


n*n  tr  Q  -  n*0  n  *  m*(D  m)  tr  (D  H)  -  m*H  m  . 

***  *v  ***  M»***  ts*  ***  ******  ***  ****** 


(3.15) 


Let  us  define  a  matrix  [L . ^ ]  in  terms  of  the  ei , 
ting 

Lii  =  XT2 (S-,  +  e2  +  s3 -  ei )  (no  sum) 

L12  *  ^2}  *  ?  ^X1  61  +  X2  s2  +  w3^  ’ 

L23  s  L32  =  \  (x22&2+  X32b3  +  w1^  ’ 

L31  s  L13  =  7  ^3^®3  +  Xl2sl  +w2^  * 


> 


in  (2.8)  by  set- 


►  (3.16) 

j 


A  considerable  amount  of  algebra,  based  on  (2.7),  (2.1),  (2.8),  and 
involving  (2.3),  (1.5),  (3.16),  enables  us  to  deduce 


m*  (D"^m)  tr(D  H)  -  m*Hm  «  2K(X;z)  ,  (3. 17) 

/■w  rv  rw  /v  r** 

provi ded 

K(i;2>  *  Luz(2d  •  2k  =  "k  •  <3-18>’ 

Consequently,  applying  (1.31),  (2.10)  to  the  first  inequality  in  (1.49), 
and  (3.15),  (3.17)  to  the  second,  we  obtain 

E(x;z)  >  0  ,  K(x;z)  >0  V  z  6 A  ,  (3.19) 

^  M  A* 

as  a  set  of  conditions,  equivalent  to  (1.49),  that  are  necessary  and 
sufficient  for  strong  ellipticity. 

We  show  presently  that  (2.40)  with  (a)  replaced  by  >  0  holds 
if  and  only  if  (3.19)  holds.  Indeed,  suppose  (3.19)  is  true.  The  first 
of  (3.19)  coincides  with  (2.15)  and  hence,  as  shown  in  Section  2,  im¬ 
plies  (2.40).  Moreover,  from  the  second  of  (3.19)  we  draw  in  particular 

K(a;?,)  >  0  , 

where  £|  ,  introduced  in  (2.13),  is  the  vertex  of  A  on  the  z^-axis 
(see  Figure  1).  Thus  (3.19)  necessitates 

L'J  1  =  (Bg  +  Bg)  >  0  » 

which,  together  with  (a)  of  (2.40),  at  once  requires 

‘Note  that  because  of  (3.16),  (2.8),  (1.4),  and  (1.14),  the  L,.  are 
solely  functions  of  the  xk  . 
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6f  >  0  .  (3.20) 

Therefore,  (3.19)  in  fact  implies  (2.40),  (3.20),  and  hence  (2.40)  with 
(a)  replaced  by  (3.20). 

Conversely,  suppose  (2.40)  and  (3.20)  hold.  Earlier  in  this  sec¬ 
tion,^  we  have  shown  that  (2.40)  implies  the  first  of  (3.19).  We  now 
substitute  for  w^  from  (2.41)  into  (3.16)  and  find  -  bearing  (1.5)  in 
mind  —  that  K(a;z)  ,  defined  in  (3.18),  admits  the  representation 

=  61 (x2iz2  "  a3  z3^  +  MA3^z3  '  V2!  ^  +  e3'Al  Z1  “  X2^z2^ 

■*■[^1  6^  ^  X2  ^  ^3  *1*2 

+  Cx22e2  +  ^32e3+X161(a)l  +  2)]  Z2Z3 

+  [x3283  + x’2B1  +  X202(o)2  + 2)]  Z3Z1  .  (3.21) 

Finally,  on  account  of  (3.20)  and  (b)  of  (2.40),  the  sum  of  the  first 
three  terms  in  (3.21)  is  positive  on  A  except  at  the  interior  point 

O 

z  defined  by 

zi  =  MA1  +  A2  +  A3^ 

whereas  the  sum  of  the  remaining  three  terms  is  positive  on  A  except  at 
the  vertices  c.,  .  Hence,  (3.20)  and  (b)  of  (2.40)  imply  the  second  of 
(3.19).  This  completes  the  proof  of  the  claim  that  (2.40),  (3.20)  hold 
if  and  only  if  (3.19)  is  true. 

^See  the  proof  of  (3.1). 
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The  foregoing  results  entitle  us  to  conclude  that  a  set  of  condi¬ 
tion^  necessary  and  sufficient  for  the  strong  ell ipti city  of  the  system 
of  partial  differential  equations  (1.10)  at  a  solution  ( jj, p )  and  a  point 
x  ,  is  given  by 


0.  >  0  , 

(a)'! 

u1  +  2  >  0  , 

O 

■> 

(b) 

[/p2(u>2  +  2) 

+  /p3(a>3+  2)  ]  , 

o 

[/p3(<i>3  +  2) 

+  /P1  ((D-|  +  2)  ]  , 

o 

►  (c) 

[/p-|  (o).|  +  2) 

+  /p2(<^2  +  2)  ] 

J 

i 

>  0.22) 


J 


For  convenience,  we  cite  here  from  (2.41),  (2.42)  that 

Pi  =  6i/xi ei s263  ’  “1  *  W1/Ai Bi  (no  sums)  , 

ei  =  ( 1^  slg)  +  A?  Wgd-j.Ig)  » 

wi  -  A?(I1  -  X^)61  -*•  2X?(  -  2I2  -  xf  -  2A~2)CW11  (I1  *I2) 

*  2a?  MjgUl  ,I2)  +  **  W22(Ii  ,I2)]  (no  sum)  , 


(3.23) 

►  (3.24) 


j 


where  A .  *x.(x)  are  the  local  principal  stretches  associated  with  u 

I  I  r*  ^ 

and  I  *  I  (a(x))  .  With  the  help  of  (3.23),  the  system  of  inequalities 
(3.22)  is  immediately  found  to  be  equivalent  to  the  more  tractable  set 


T» 

>li 

V 

O 

«• 

(a) 

wi  +  ^^i0i  >  0  (no  sum)  ’ 

(b) 
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x‘2(Wl  -  2A1e1)  <  [/a"2(w2  +  2x2B2)  +  w3  +  2x383)]2  , 
\~2Z(w2- 2\2B2)  <  [,/a“2(w3  + 2a383)  +  +  2x1b1)  ]2  ,>(c) 

X"2(w3- 2X383)  <  [Jx"2(w]  +2X181)  +  yx"2(w2  +  2X282)  ]2  . 


(3.25) 


Here  are  defined  in  terms  of  the  derivatives  of  W  by  (3.24) 

★ 

or  alternatively  in  terms  of  the  derivatives  of  W  through  (2.53), 
(2.54). 

In  order  to  bring  out  the  invariance  of  (3.25)  under  permutations 
of  (x1,x2,X3),  we  define  functions  xi  by 


9  >^2^  =  W^1  ’^2^  +  ^*1  *2  ® ^*1 f *2^  * 


x3^1  =  *3  ^^3® ^ * ^2 ^  ~  ^ (^1  * ^2)] 

+  +  J^2  x2^3’*l  ^ ^  *  ^3  =  (A-|A2) 


(3.26) 


where  e  and  w  are  given  by  (2.54);  thus 


s(  A-j  ,x2)  - 


i  ★  2 

g  w11(x,x“t)  for  A-j  =  A2  =  A 


r\  x2(x2  -  x2)"1W1(x2,x3)  for  X1J(X2,  x3=(x1x2)'1, 


/ 4  A  W-ji  (x,x  )  for  X1  =  x2=  x 

W  (  X 1  9X3)  =V 

\1(x^-x2)-1w1(x2,x3)  +  ^  x^2w11(x2,x3) 


for  x ^  f  Xg  ,  X3 —  ( x^  A2 )  . 


(3.27) 


In  view  of  (2.53),  the  set  of  inequalities  (3.25)  may  now  be  written  as 
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x.j  ,^2)  >  0  ,  ^i^2’^3^  >  0  ,  xi(x3,x1)  >  0  .  (3.28) 

Moreover,  from  (3.26)  and  the  first  two  of  (2.59),  we  infer 

x^x-j.Ag)  =  x^ (AgfX^ )  .  (3.29) 

By  virtue  of  (3.29),  the  set  of  strong  ellipticity  conditions  (3.28)  - 
like  the  system  of  ordinary  ellipticity  conditions  (2.58)  —  is  fully 
symmetric  with  respect  to  the  principal  stretches. 

We  now  prove  that  the  conditions  for  ordinary  and  strong  el  1 ipticity 
coalesce  for  an  isotropic  material  of  the  type  under  consideration  that 
obeys  the  Baker-Ericksen  inequalities  (1.19)  and,  in  addition,  has  a^ 
positive  shear  modul us  at  infinitesimal  deformations. 

Clearly,  in  the  case  of  distinct  A i  ,  the  inequalities  (1.19)  - 
in  view  of  (3.24),  (1.5)  -  imply  (3.20)  and  hence  the  equivalence  of 
(2.40)  and  (3.22).  On  the  other  hand,  suppose  two  of  the  principal 
stretches  coincide  but  are  different  from  the  third,  say  ^  \y 
Then  (1.19),  (3.24),  (1.5)  give  =  62  >  0  »  and  ^gain  (2.40)  holds  if 
and  only  if  (3.22)  is  true. 

Finally,  suppose  all  three  A.  are  equal,  so  that  A^  =  l.  In 
this  degenerate  instance,  (2.53),  (2.54),  and  (3.23)  yield 

6rrwi =  s"n(1’1}  *  «i-2  •  0.30) 

Thus,  conditions  (2.40)  and  (3.22)  at  present  respectively  reduce  to 

W^d.l)  t  0  and  W^d.l)  >  0  .  (3.31) 

One  easily  confirms  that  if  u  is  the  shear  modulus  of  the  material  at 
infinitesimal  deformations. 
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u  =  2[W1(3,3)  +  W2(3,3)]  =  |  W^O.l)  . 


(3.32) 


Accordingly,  when  p  >  0  ,  the  two  conditions  in  (3.31)  are  equivalent. 
This  concludes  the  proof  of  the  preceding  claim  concerning  the  equiva¬ 
lence  of  ordinary  and  strong  ellipticity. 

It  is  of  interest  to  note  upon  substitution  from  (3.24),  (1.4), 
(1.5)  into  (a),(b)  of  (3.25),  that  these  inequalities  hold  if  and  only 
if 


W-|  +  A?  W2  >  0 


^1  +  Xi  ^2  +  2(  I.j  -  A?  -  2x  •"*  )(W-ji  +  2x?  W^2  +  A4  W22)  >0  (no  sum) ,  J 


(3.33) 


and  consequently  one  recovers  the  necessary  conditions  for  strong  ellip- 
ticity  cited  by  Sawyers  [12]. 1  To  see  that  the  additional  restrictions 
(c)  of  (3.25)  are  in  general  independent  of  (a)  and  (b),  consider  the 
elastic  potential  defined  by 


W(I1 ,12)  =  !i  +TF  I2  -  l\  ’  W  (3>34) 

and  a  deformation  with  local  principal  stretches 

X1  =  \  ’  X2  =  2  ’  x3  =  1  -  (3.35) 


In  this  instance,  (a)  and  (b)  of  (3.25)  are  satisfied,  while  the  first 
inequality  in  (c)  is  violated. 


T 


See  inequalities  (4.1)  of  [12]. 


) 
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4.  El  1 i pt i city  conditions  for  special  classes  of  deformations  and 
materials. 

In  applying  the  ellipticity  conditions  established  in  the  preceding 
sections  to  particular  classes  of  deformations  and  hyperelastic  solids, 
we  shall  take  for  granted  that  the  material  obeys  the  Baker-Ericksen 
inequalities  (1.19)  and  has  a  positive  shear  modulus  at  infinitesimal 
deformations.  Hence,  there  will  be  no  need  to  distinguish  between  the 
criteria  for  ordinary  and  those  for  strong  ellipticity. 

We  consider  first  a  deformation  with  local  principal  stretches 


If  t.  is  the  principal  true  stress  correspond!' ng  to  ,  induced  by 
this  deformation,  (1.17)  gives 

t-|  =  t2  .  (4.2) 

Conversely,  (4.2),  (1.17),  and  (1.4),  (1.5)  -  because  of  (1.19)  —  imply 
(4.1).  Thus,  the  special  deformation  at  hand  is  associated  with  a 
locally  axlsymmetric  state  of  true  stress.  Further,  by  virtue  of  the 
hydrostatic  pressure  field  p  ,  this  state  may  in  particular  be  one  of 
local  uniaxial  stress,  in  which  case  tiSt2=0;  alternatively,  it  may 
be  a  state  of  local  equi-biaxial  stress  with  t3  =  0. 

From  (4.1),  (2.53),  (2.54),  and  the  first  two  of  (2.59)  follow 

-2  -2 

81='e2sS(A  .A)  ,  w1=w2  =  w(x  ,a)  ,  2a383  =  w3  =  w(a,a)  .  (4.3) 

With  the  aid  of  (4.3),  the  ellipticity  conditions  (3.25)  for  a 
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deformation  locally  characterized  by  (4.1)  are  readily  found  to  reduce 
to 

S]  >  0  ,  >  0  ,  w^+2A8-|  >  0  .  (4 . 4 ) 1 

We  show  now  that  the  first  of  these  inequalities  is  automatically  satis¬ 
fied  under  our  present  assumptions.  Indeed,  if  A f  1  ,  the  Baker- 
Ericksen  inequalities  (1.19)  and  the  first  of  the  relations  (3.24)  at 
once  give  8-j  >  0  ;  if  A  =  1  ,  in  turn,  then  B-j  >  0  follows  from  (3.30) 
and  (3.32)  with  u >  0  .  Thus,  bearing  in  mind  (3.24),  (1.4),  and  (4.1), 
one  arrives  at 


Wi  +  A  4W2  >  0  , 

X4(W1  +  A2W2)  +  2(a3- 1)2(W1] +2a2W12  + A4W22)  >0,  J 


(4.5) 


with 


'VVW  •  “a6  =  “ae(IrV  *  V2*2-4 


,  I2  =  2a 


'2  +  a4, 


(4.6) 


as  necessary  and  sufficient  eilipticity  conditions  in  the  circumstances 
under  consideration. 

It  is  also  useful  to  cast  (4.5),  or  alternatively  the  last  two 
inequalities  in  (4.4),  in  terms  of  the  partial  derivatives  of  W  . 

We  note  from  the  first  relation  in  (a)  of  (2.49)  and  the  second  in  (b) 


that 


W'j  (a, a”2)  =  0  V  a  >  0  . 


(4.7) 


Substituting  from  (3.24),  (4.1)  into  (4.4)  and  invoking  (1.4),  one 
recovers  the  inequalities  (4.8)  of  Sawyers  [12],  who  observes  that  the 
latter  are  both  necessary  and  sufficient  for  strong  eilipticity  so  long 
as  the  deformation  locally  conforms  to  (4.1). 
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Thus  and  by  (4.3),  (4.1),  (3.27),  the  last  two  of  (4.4)  are  equivalent 
to 

W^U.A-2)  >  0  ,  2W1  (X,x)  +  A(1  +  X3)  W]  1  ( A , A )  >  0  ,  (4.8) 

so  that  these  two  inequalities  are  also  necessary  and  sufficient  ellip- 
ticity  criteria  when  the  local  principal  stretches  satisfy  (4.1). 

Consider  next  a  deformation  corresponding  to  a  state  of  local  plane 
strain  with  principal  stretches 

=  X  ,  -  X  ,  X ^  —  1  •  (4.9) 

In  this  instance  our  assumptions  (1.19)  and  y>0  are  easily  seen  to 
imply  e.j  >  0  .  In  fact,  if  Xj*  1,(1. 19)  and  (3.24)  immediately  yield 
6- >  0  ,  whereas  for  A  =  1  the  Bi  are  positive  as  a  consequence  of 
(3.30)  and  (3.32).  Thus,  for  local  plane  strain,  the  only  ellipticity 
conditions  to  survive  are  (b)  and  (c)  in  (3.25). 

Abeyaratne  [6]  deduced  necessary  and  sufficient  conditions  of  or¬ 
dinary  ellipticity  for  a  deformation  of  global  plane  strain  with  a 
displacement  field  obeying 

ua,3  =  0  -  u3  =  0  '  (4-10) 

Such  a  deformation  gives  rise  to  (4.9)  throughout  the  body.  The  condi¬ 
tions  of  ordinary  ellipticity  arrived  at  in  [6]  are 

W'(I)  f  0  ,  1  +2(1-2)  [W"(I)/W’(I)]  >  0  ,  (4. II)1 

^See  inequalities  (3.21)  in  [6]. 
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where 

I  =  - 1  =  I2-  1  =  x2  +  a-2  ,  W(  I)  =  W(  I  1  ,  1  +  1)  .  (4.12) 

We  now  verify  that  (4.11)  is  at  present  implied  by  (b)  and  (c)  in  (3.25), 
as  must  be  the  case  when  the  Baker-Ericksen  inequalities  (1.19)  are 
assumed  to  hold  and  u>  0.  Indeed,  one  gathers  from  (3.24),  (4.12), 
and  (4.9)  that 

W'  ( I )  =  83  ,  W(I)  +  2(I- 2)  W"(I)  =  A2(x2+ir2(w3+2x3e3).  (4.13) 

But  w^ + 2*363 >0  according  to  (b)  in  (3.25),  while  63 >  0  under  the 
assumptions  just  mentioned,  whence  (4.11)  follows.^ 

The  necessary  and  sufficient  ellipticity  conditions  appropriate  to 
a  local  state  of  plane  strain  governed  by  (4.9)  are  in  general  far  more 
restrictive  than  (4.11)  since  the  former -involve  all  three  inequalities 
in  (b)  of  (3.25),  as  well  as  the  three  inequalities  in  (c).  These  addi¬ 
tional  restrictions  stem  from  the  fact  that  the  local  hypothesis  (4.9) 
admits  a  larger  class  of  deformations  than  does  the  global  requirement 
(4.10);  furthermore,  (4.11)  —  with  (1.19)  in  force  —  preclude  potential 
discontinuities  in  the  relevant  normal  derivatives  of  the  displacements 
and  pressure  merely  across  cylindrical  surfaces  with  generators  parallel 
to  the  principal  direction  of  the  deformation  tensor  C  that  is  asso¬ 
ciated  with  X3  =  1  . 

^It  should  be  noted  that  (4.11),  as  expected,  is  also  implied  by  the 
conditions  of  ordinary  ellipticity  (2.40)  even  if  the  Baker-Ericksen 
inequalities  fail  to  hold  or  if  u  <.  0  . 
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We  proceed  now  to  special  classes  of  materials  and  in  this  connec¬ 
tion  consider  first  elastic  potentials  ’at  depend  on  a  single  deforma¬ 
tion  invariant.  Suppose 


W(IrI2)  =  «(!,) 


(4.14) 


We  show  presently  that  in  this  event,  (a)  and  (b)  imply  (c)  in  the 
ellipticity  conditions  (3.25)J  From  (4.14)  and  (1.4),  (1.5),  (3.24) 
we  draw 


=  W'Uj)  , 

xT2(Wi  t  2xi6i )- (I, -*f±  2a:1  )CW'(I1)  +  2(I1  -xf  ?  2X-1  JWd, )] 

(no  sum)  . 


V  (4.15) 


Assume  now  that  (a)  and  (b)  in  (3.25)  hold  true,  or  equivalently 


wd,)  >  o  , 


W(I1)  +  2(I1  -  xf  -  2X^])  W"(I1)  >  0  . 


(4.16) 


With  a  view  towards  inferring  the  first  inequality  in  (c)  of  (3.25),  we 
note  that  the  latter  can  be  written  as 


S  <  T1/2  ,  (4.17) 


provided 

S  =  A^2(w1-2xTe1)- a'2(w2  +  2a262)  -  a’2(w3  +  2a303)  ,  j 

>  (4.18) 

T  «  4a2(w2  +  2a2B2)  (w3  +  2a3s3)  .  J 


Clearly,  T  >  0  under  our  current  hypotheses. 


^One  may  show  similarly  that  (a)  and  (b)  imply  (c)  in  the  ordinary 
ellipticity  conditions  (2.40)  whenever  (4.14)  holds. 


On  account  of  (4.18)  and  (4.15),  (1.4),  (1.5), 


S  =  -2(x1  +  x2)(x]  +  x3)[W'(I1)  +  2(x1  -  x2)(x1  -  x3)  W,,(I1 )]  , 
T  =  S2  +  U  , 


Y  (4.19) 


where 

U^fx^Xj)2  (X1+X3)2  (X2-X3)2W'(I1)W,,(I1)  .  (4.20) 

In  order  to  establish  (4.17),  it  is  convenient  to  observe  that 


either  x2=x3  or  x^x^,  W"(I^)<.0  or  x2  f  x3,  WM( I-|  )  >  0.  (4.21) 

If  the  first  of  these  three  eventualities  holds,  (4.16),  (4.19)  imply 

S=-2(x]  +x2)2[W'(I1)  +  2(x1  -x2)2  Wd,)]  <  0  , 

and  hence  (4.17)  follows.  On  the  other  hand,  the  second  alternative  in 
(4.21)  justifies 

2(x1  -  x2)(x1  -  x3)W,,(I1)  >  min[2(x1  -  x2)2W,,(I1),  2(x]  -  x3)2W"(I1)]  , 

so  that  (4.16),  (4.19)  again  yield  S<0  and  thus  (4.17).  Finally, 
the  last  alternative  in  (4.21),  together  with  (4.20),  (4.19),  and  the 
first  of  (4.16),  gives  U>0  and  T  >S  ,  which  confirms  (4.17)  in 
this  case  also.  Strictly  parallel  arguments  enable  one  to  infer  the 
remaining  two  inequalities  in  (c)  of  (3.25). 

As  is  now  clear,  (4.16)  constitutes  necessary  and  sufficient  ellip- 
ticity  conditions  when  the  elastic  potential  obeys  (4.14).  Moreover, 
the  Baker-Ericksen  inequalities  and  y  >  0  imply  the  first  of  the  two 
inequalities  (4.16.) 
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If,  instead  of  (4.14), 

W(IrI2)  =  W(I2)  ,  (4.22) 

one  finds  similarly  that  (a ) , (b )  once  again  imply  (c)  in  (3.25)  and  is 
led  to 

W(I2)>0,  W,(I2)  +  2(I2-a:2-2x.)  W"(I2)  >  0  ,  (4.23) 

as  necessary  and  sufficient  ellipticity  conditions. 

It  should  be  mentioned  that  the  ellipticity  conditions  (4.16)  and 
(4.23),  appropriate  to  (4.14)  and  (4.22),  respectively,  recover  earlier 
results  reported  by  Sawyers  [12].1 

Consider  now  the  special  case  of  a  Mooney-Rivl in  material,  for 

which 

W(IrI2)  =  C]  (I1  -  3)  +  c2(I2-3)  .  (4.24) 

The  Baker-Ericksen  inequalities  (1.19)  necessitate 

co>.0  ,  p  =  yc2  >  0  ,  (4.25) 

and  with  (4.25)  in  force,  the  ellipticity  conditions  (3.25)  are  found 
to  be  automatically  satisfied  for  all  principal  stretches.  Thus,  for  a 
Mooney-Rivl in  material  that  obeys  the  Baker-Ericksen  inequalities, 
ellipticity  obtains  at  all  deformations. 

Finally,  as  an  illustrative  example,  we  apply  the  ellipticity  con¬ 
ditions  established  in  this  paper  to  a  class  of  hypothetical  materials 

^See  Section  5  of  [12]. 


I 


introduced  by  Ogden  [13].  The  associated  mechanical  response  is  gov¬ 
erned  by  an  elastic  potential  of  the  form 

W(ArA2)=|  (x“+a“+a“- 3)  ,  A3=(A1A2)"1,  c  >  0  ,  a  >  0  .  (4.26) 

According  to  (3.27)  and  (4.7)  one  has  at  present 


6  (  A-j  » A2  )  = 


xa  ^  for  =  x2  =  A 
f  (x^-x“)/(x^-x2)  for  X-|  f  x2  , 


(A-jAg)  [w(X^  ,X2)  +  2(A^X2)  6(A.|,A2)] 

=  f  A“(l+6)-1[(a.l)(l  +  6a+1)+(c+l)(5  +  6a)]  , 


(  A 1  A2  )  [w(  A-j  ,  A2  )  -  2  (  Xi  A2  )  B(Ai  ,  A2  )  ] 


0  for  x^  =  x2 


(4.27) 


‘  '  2  l  (a  - 1  J(<s“+ 1  )J 


for  Ai  f  a2  , 


where  6  =  x-|/x2.  With  the  aid  of  (4.27)  and  by  recourse  to  (2.53),  one 


can  show  that 


B.  >  0  , 


A^c(wi  +  2xis1 )  >  |  (a“  +  a2  +  A3  -  A?)  (a  -  1 )  (no  sum) 


A^t(wi  -  2a.b1)  <  |  (x“  + x“+ x“- a?)  max[0,  a- 1]  (no  sum)  . 


(4.28) 


The  first  of  (4.28)  implies  the  Baker-Ericksen  inequalities  (1.19),  and 


i 
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(4.26) -together  with  (3.32)  — gives  u  =  ca/2  >  0  ;  further,  (a)  in  (3.25) 
is  satisfied  for  all  x^  and  all  a>0  .  Moreover,  when  a>.l  ,  we 
easily  see  that  (b)  and  (c)  in  (3.25)  also  hold  for  all  x^  .  In  con¬ 
trast,  as  may  be  shown  by  means  of  the  second  of  (4.27)  and  the  last 
of  (4.28),  for  0 <  a  <  1  ,  (b)  and  (c)  hold  if  and  only  if 

e  <  X^/Xj  <  e  ^  »  (4.29) 

provided  e  =  e(a)  is  the  unique  solution  of  the  equation 

(1 +a)(e  +  £a)  =  (1  -  a)(l +ea+1)  ,  0<e<l  .  (4.30) 

Therefore,  for  an  elastic  potential  obeying  (4.26)  with  ail,  el  1  i pti - 
city  prevails  at  all  deformations,  whereas  if  0<a<l  ,  the  inequali¬ 
ties  (4.29)  are  necessary  and  sufficient  for  ellipticity.1 

The  inequalities  (4.29),  which  are  contingent  upon  0<a<l  , 

describe  the  interior  n  of  a  convex  hexagonal  pyramid  in  the  space  of 

2 

the  principal  stretches,  the  origin  of  which  coincides  with  the  vertex 
of  this  pyramid.  The  intersection  z  of  Q  with  the  "isochoric  surface" 
*1A2A3*1  represents  the  ellipticity  domain  in  this  instance. 

Evidently,  (4.30)  can  be  solved  explicitly  if  a  =  j.  In  this 
special  case,  the  ellipticity  conditions  (4.29),  upon  elimination  of 
A j  ,  may  be  written  as 

e<X1X21<e"1,  e  <  X-j  *2  <  E"^  *  e<^x2<e"1,  £  =  7-4/3^0.0718. 

(4.31) 

^ Ogden  [13],  on  the  basis  of  tests  performed  on  vulcanized  natural  rub¬ 
ber,  suggests  the  choice  a  =1.2  and  c*6.8  kg/cm2. 

2 

See  the  interpretation  of  inequalities  (3.1)  of  [1]. 
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Conditions  (4.31)  define  the  open  region  z^2  in  the  (x-j  ,x2)-plane 
depicted  as  the  shaded  area  in  Figure  2.  This  region  is  the  projection 
onto  the  (x^ ,x2)-plane  of  the  ellipticity  domain  Z  ;  the  projections 
z23  and  z31  of  1  onto  the  other  two  principal -stretch  planes  are  once 
again  represented  by  the  shaded  region  in  Figure  2,  provided  (x^,x2) 
in  this  diagram  is  replaced  by  (x^.x^)  or  by  (x^.x^  . 

It  is  of  interest  to  examine  the  response,  to  certain  basic 
homogeneous  deformations,  of  the  material  characterized  by  the  elastic 
potential  (4.26).  For  a  pure  homogeneous,  volume-preserving  deforma¬ 


tion  of  the  form 


yi  =  x^  (no  sum)  ,  x-jX2x3  =  1 


(4.32) 


one  gathers  from  (1.7),  (1.16),  the  first  of  (2.46),  and  (4.26)  that 


'ij  ’  3  0  • 


while 


xiai13Tii  *  tbb"t33=  abVW  “  c^xb'AP  ^no  sum  over  1  or 

(4.33) 

In  particular,  for  a  pure  homogeneous  deformation  of  (i)  uniaxial 
stress  or  (ii)  equi-biaxial  stress  or  (iii)  plane-strain  uniaxial 
stress  ("pure  shear"),  (4.33)  yields  respectively: 


(0  Tn  *  t22 


=  0  ,  Xo33(x)  =>  t33(x)  =  c(xa  -  x_a/2),  X,  =  X,  X-,  =  X,  =  X-1/2; 


V1  2 


(ii)  t33*0,  Xa22(x)  =  t22(x)  =  c(xa-  x"2a)  ,  X1  =  X2=  X  ,  X3=  x'2  ; 


(iii)  t22  =  0  ,  xa11(x)*r11(x)  =  c(xa-x‘0t)  ,  X]  =  X"1  =  X  ,  x3=  1  . 


(4.34) 
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Next,  consider  a  deformation  of  simple  shear  given  by 

y^  ~  x-j  +  kx2  >  y 2 =  x £  ,  y3  =  x3  >  (4.35) 

where  k  is  the  amount  of  shear.  The  corresponding  principal  stretches 
obey 

■  x'1  =  \  (k+^  +  4),  x3=l  ,  k=x1-x'1  ,  (4.36) 

and  the  in-plane  true  shear  stress  here  is  found  to  satisfy 

x12  =  t(x)  =  c(xa- x“a)  /  (x  +  x-1)  ,  x=x]  .  (4.37) 

When  a  >.1 ,  the  normal  stresses  o33(x),  t33(x)  and  a22(x),r22(x), 
as  well  as  (x),t,j(x),  arising  in  the  pure  homogeneous  deformations 
(4.34),  are  monotone  increasing  functions  of  X  for  0  <x<°°.  When 
0<a<l  ,  on  the  other  hand,  this  monotonicity  property  is  retained  by 
the  true  stresses  T33(x),  t22(x),  t-|]  M  »  but  the  associated  nominal 
stresses  °33(x),  a22(x),  a^-j(x)  are  steadily  increasing  merely  for 
0<  x  <  X*  ,  possess  a  maximum  at  x  =  x*  ,  and  monotonically  decrease  to 
zero  on  x*<x<»  ,  where  from  (4.34), 

(D  **  -  [(2  +  c)/2(1-o)]2/3°  >  l  , 

(ii)  **  *  [(l  +  2a)/  (l-a)]1/3a  >  1  , 

(iii)  **  *  C(1  +  cc)  /  (1  -  a)]1/2ot  >  1  , 

respectively.  The  graphs  of  a33(x)  and  t33(x)  appropriate  to  case  (i) 
of  uniaxial  stress  are  shown ^  in  Figure  3  for  the  three  values  ct=  3/2  , 

^In  the  interest  of  clarity,  the  curves  for  t«(x)  have  been  omitted 
for  0  <  x  <  1  . 


>  (4.38) 
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a=l  ,  a  =  1/2  .  The  graphs  of  c22(x),  t22(x)  and  corre¬ 

sponding  to  the  cases  (ii)and  (iii),  are  qualitatively  quite  similar. 

As  far  as  the  simple-shear  deformation  (4.35)  is  concerned,  (4.37) 
reveals  t(a)  to  be  steadily  increasing  over  the  entire  range  0<x<®, 
provided  a>.l  .  When  0<a<l  ,  however,  t(x)  increases  monotoni- 
cally  only  for  e^2<a<e~^2,  with  e  =  e(ct)  given  by  the  root  of 

(4.30);  for  this  range  of  the  parameter  a  ,  t(x)  has  a  minimum  at 
1/2  -1/2 

X  =  e  '  ,  a  maximum  at  X  =  e  '  ,  and  is  strictly  decreasing  for 

0<X<E1/2  and  e"^2<X<®  .  The  graphs  of  t(a)  corresponding  to 
a  *  3/2  ,  a  =  1  ,  and  a=  1/2  are  displayed  in  Figure  4. 

It  is  evident  from  the  foregoing  discussion  of  the  three  pure 
homogeneous  deformations  introduced  in  (4.34)  that  when  0<a<l  ,  the 
nominal  normal  stresses  singled  out  there  are  no  longer  invertible  func¬ 
tions  of  the  appropriate  principal  stretch  X  over  the  range  1  <  x <  , 

which  corresponds  to  an  elongation.  No  such  loss  of  invertibility  is 
encountered  in  contraction,  i.e.,  for  0  <  x  < 1  .  Furthermore,  el  1 i pti - 
city  is  lost  when  0<a<l  in  both  extension  and  contraction,  and  as 
illustrated  by  Figure  3  -  which  pertains  to  uniaxial  stress  —  the  maxi¬ 
mum  of  033(A)  occurs  in  the  interior  of  the  range  of  ellipticity.1 

As  regards  the  behavior  in  simple  shear  of  the  special  material 
under  discussion,  we  note  that  here,  for  0<a<l  ,  the  range  of 

^n  case  (ii),  which  refers  to  equi -biaxial  stress,  invertibility  of 
o2?(x)  is  lost  once  again  only  in  extension,  but  this  loss  occurs  out¬ 
side  the  range  of  ellipticity  for  0<  a  <0**0.714  and  falls  inside 
this  range  for  a* <  0  <  1  . 


ellipticity 


[e(a)]1/2  <  X  <  [e(a)]-1/2 

precisely  coincides  with  the  range  of  X  over  which  the  shear-response 

function  t(x)  is  monotone  increasing,  and  the  two  extrema  of  t(x) 

occur  at  the  endpoints  of  this  range  (see  Figure  4). 

It  is  now  also  apparent  that  for  the  material  characterized  by 

(4.26),  ellipticity  prevails  at  a  local  state  of  plane  strain  with 
“1  ° 

x]  =  x2  =  X,  x^*!  if  and  only  if  the  shear-response  function  r(x)  has 
a  positive  slope  at  x  =  x  .  This  conclusion  is  in  complete  accord  with 
a  result  previously  arrived  at  by  Abeyaratne  [6].  A  related  conclusion 
had  been  reached  still  earlier  by  Knowles  [4],  who  found  that  the  dis¬ 
placement  equations  of  equilibrium  appropriate  to  global  anti-plane 
shear  of  a  class  of  incompressible  hyperelastic  materials  are  locally 
elliptic  at  a  solution  if  and  only  if  the  response  curve  for  simple 
shear  has  a  positive  slope  at  an  amount  of  shear  equal  to  the  magnitude 
of  the  local  displacement  gradient. 

The  preceding  observations  suggest  that  a  loss  of  ellipticity  in 
finite  elastostatics  always  entails  some  loss  of  invertibility  of  the 
underlying  stress-deformation  relations.  A  physical  interpretation 
along  these  lines  of  the  three-dimensional  ellipticity  conditions  estab¬ 
lished  in  this  paper  has  so  far  eluded  our  efforts. 
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FIGURE  3.  BEHAVIOR  OF  THE  SPECIAL  MATERIAL 
UNOER  UNIAXIAL  STRESS. 
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FIGURE  4.  RESPONSE  OF  THE  SPECIAL  MATERIAL 
TO  SIMPLE  SHEAR. 
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